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I. Introduction. 


In the last few years, a number of algorithms using spectral collocation methods have been 
successfully implemented to solve the incompressible Navier-Stokes equations in a domain 0 of 
or 

(1.1) - V Au + grad p + (u,V)u = f in Q , 

div u = 0 in Q , 

provided with appropriate boundary conditions : periodic. Dirichlet or mixed ones, We refer to 
[VGH] and to [CHQZ] for a detailed bibliography. In most of them, the collocation points in the 
nonperiodic directions are the nodes of a Gauss-Lobatto quadrature formula associated with the 
Chebyshev polynomials. Indeed the use of the Fast Fourier Transform allows a less expensive 
computation of the derivatives and the nonlinear terms. However, as far as we know, the only 
theoretical justifications of some of these algorithms are achieved with the Chebyshev nodes 
replaced by the Legendre ones (see e.g, , [BMM 1 ][BMM2]), The aim of this paper is the numerical 
analysis of a collocation method involving the Chebyshev nodes, for Dirichlet boundary conditions 

r\ 

and when the domain Q is the square ]-l ,1[ ; corresponding numerical experiments can be found 
in [Mo][Me]. 

The analysis we present here has already been achieved for the Legendre collocation nodes, 
in [BMM2], The extension to the Chebyshev methods presents essentially two difficulties : 

1 ) In the variational formulation of both the continuous and the discrete problems, the classical 
Sobolev spaces have to be replaced by weighted Sobolev spaces, the weight of which is the 
Chebyshev one. Several trace theorems in these spaces, as well as some regularity results for the 
Dirichlet problem for the laplacian in a square, are needed by our study ; these are recent (see 
[BM]). 

2) Due to the Chebyshev weight, the analysis of the Stokes problem - which one obtains by 
neglecting the nonlinear terms in (1,1) - and of its approximation involves a variational 
formulation of the form 

(1.2) a(u,v) + b^(v.p) = <f,v> , 

b2(u.q) = 0 

In this system, the bilinear forms b, and b 2 are distinct, in opposition to the classical 
saddle-point problem first studied in [Br], We shall consider such a problem in an abstract 
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framework and prove that it is well-posed if and only if several inf-sup conditions are satisfied : 
two for the form a and one for each of the forms b, and bg . We shall also analyze its discretization 
and state general error estimates which can be applied to a number of numerical methods. Our 
hope is that this abstract variational formulation will be used for other equations involving 
weighted spaces. 

In the particular case of the discretization of the Stokes problem by spectral methods, in 
order to obtain a well-posed problem , we have to exhibit the spurious modes of the pressure, i.e. 
the modes which cancel the discrete gradient. Thus, we derive an appropriate choice for the space 
of discrete pressures and we can prove that the inf-sup conditions are satisfied. The results we 
obtain are very similar to the corresponding ones for the Legendre methods (see [BMM2]). The 
convergence estimates can be extended to the Navier-Stokes equations without difficulty. 

An outline of the paper is as follows. Seotion 1 1 is devoted to an abstract variational problem 
with inf-sup conditions. In Section III , we state a variational formulation of the Stokes problem in 
weighted Sobolev spaces, the weight of which is the Chebyshev one. In Section IV, we study a 
Galerkin spectral method to discretize this problem. The spectral collocation method for the Stokes 
problem is thoroughly analyzed in Section V. The results are extended in Section VI to the full 
Navier-Stokes equations thanks to a fixed point theorem. 

Notations : The norm of any Banach space E is denoted by IMI^ . while <.,.> is the duality pairing 
between E and its dual spaoe E'. For any pair (E.F) of Banaoh spaces, L(E,F) represents the space of 
continuous linear mappings from E into F. We mean by A 0 B the tensorial product of any sets A 
and B in a Banaoh space, while A®^ is the tensorial product of A by itself. 

For any domain A in 1R‘* and for any real number s, we use the classical hilbertian Sobolev 
spaces H®( A) , the norm of which is denoted by ||.||j . 

In all that follows, c, c' ... are generic constants, independent of the discretization. 
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II. An abstract variational sustem. and its approximation. 


11.1. The continuous case. 

Let Xj and M, (i = 1,2) be real reflexive Banach spaces. We assume we are given three 
continuous bilinear forms, a : Xg x X, -»■ [R and bj : Xj x Mj ->■ IR (i = 1 , 2). For any given f in X', 
and g in Mj , we consider the following problem -.Find (u,p)/n Xg x such that 


( 11 . 1 ) 


Vv€X, , a(u,v) + b,(v,p) 3 <f,v> 
V q € Mg , bg(u,q) = <g,q> . 


In order to study this problem, let us introduce the linear operators A eLCXg ,X^) and 
Bj €/.(X| ,Mj') (i = 1 , 2) assooiated with the forms a and bj by the relations 

(11.2) V u € Xg , V V € X^ , <Au,v> = a(u,v) , 

(11.3) Vu€Xj,Vq€Mj, <BjU,q> = bj(u,q) ; 

we denote by Bj'' €L{t^■ ,X:) (i = 1 , 2) the adjoint operator of Bj . 

For any g in Mj' (i = 1 , 2) , we define the closed affine space 

(11.4) Kj(g) = { V € Xj ; V q € Mj , bi(v,q) = <g,q> } , 

and we note that the subspace Kj = Kj(0) is the kernel of the operator Bj . Moreover we introduce 
the linear mapping TT : X!, -> K, defined for each f e X^ by 

(11.5) V V 6 K, , <rtf,v> = <f,v> 

Finally, we denote by K° the polar set of K, i.e., thekernel of the mapping TT. 

Let us define the operator A : Xg x M, ^ X:, x Mg by the relation 

(11.6) A(u,p) = (Au + b|p , Bg u) ; 
then (11.1) is equivalent to 

(11.7) A(u,p) = (f,g) , 


We can prove the following theorem : 

Theorem 11.1 : The operator A is an isomorphism from Xgx M, onto X', x Mg if and only if 
the following conditions (CQ)snd (Cj) (i = 1 , 2) are satisfied 
(Cq) TTA: Kg K'j is an isomorphism , 

(Cj) there exists a constant gj > 0 ,such that for any q € Mj , llBj'^qllj^; > ^jllqll^. • 

Proof : We assume first that conditions (Cq) and (Cj) (i = 1 , 2) are satisfied, and we prove that A 
is an isomorphism. Given f in X', and g in Mj , let us observe that by (Cg) and [Br, Thm 0.1 ] there 
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exists u in Xg such that BjU = g and Hul^ ^ by (Cq) there exists a unique Uq 

in <2 such that TTAUq = TTf - TTAu , and one has 

IKIIx2<c( |lf||x; + ilg|lM2 ) ■ 

Let us set u = U + Uq . The element f - Au belongs to kJ . hence by (C,) and the Closed Range 
Theorem, there exists a unique p in M, such that Bjp = f - Au . Thus we have proved that A is 
onto. It is straightforward to check that A is one to one, and since it is continuous because so are 
A, b| and Bj , we conclude that A is an isomorphism by using the Open Mapping Theorem, 

Conversely, let us suppose that A is an isomorphism, For each q in M, , we have 
A(0,q) = (B|q,0); hence, from the continuity of A’’ , we obtain 
«q||„, <l|A-'UByq|lx; , 

which is (C,) . On the other hand, for each g in Mj . let us set (w.q) = A''(0,g), Then the 
mapping : g->w is continuous from M^ into X, and is such that B 2 W = g. Thus we have (C 2 ) , using 
[Br, Thm 0,1], Now we prove that TT A is an isomorphism. Given f in , let us set 
(w,q) = A"’(f,0). Then w belongs to Kj and satisfies 
HAw = nf- nB|q = nf = f , 

since nB| is equal to 0. Thus TIA is onto. Finally, let w in Kg be such that TTAw = 0; by (C,) and 
the Closed Range Theorem, there exists q in M, such that B|q = - Aw. Thus A(w,q) is equal to 
(0,0) , hence, w is equal to 0, i.e,, TTA is one to one. This concludes the proof of the theorem. 

Remark 11.1 : By well-known results of functional analysis, we can express condition (Cq) in a 
variational form. Precisely, (Cq) is equivalent to the following condition (see, e.g. [Ba]) : there 
exists a constant o(^ > 0 such that 

(11.8) Vu€K 2 , sup 

V € K, ll''llx, 2 

and 

(11.9) Vv€K,\{ 0} , sup a(u,v) > 0 . 

u € Kg 

By the Open Mapping Theorem this is also equivalent to the existence of a constant oCg > 0 such 
that 

(11.10) Vv€K, , sup ^a 2 llv|lx, 

u € Kg Il'^liXg ' 

and 

(11.11) VueKpXlO) , sup a(u,v) > 0 , 

V € K, 
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If K, (or K 2 ) are finite dimensional spaces, then the relation (11.9) or (II, 1 1 ) can be replaced by 
the requirement that 
( 11 . 1 2 ) dim K, = dim . 

Similarly, we can write the condition (Cj) (i = 1 , 2) equivalently as follows : there exists a 
constant > 0 such that 

vTx, ■ ■ 

These are the forms under which the conditions (Cq) and (Cj) (i = 1 , 2) are usually checked in the 
applications. 


Following the proof of Theorem 11.1 one can easily estimate the norm of the inverse 
isomorphism A"’ in term of the constants associated with the forms a and bj . 

Denoting by if the norm of a ; 


a(u,v) 


(11.14) X = sup 

U € X2 , V € X ^ 

we have the following result. 

Corollaru 11.1 : Assume that huootheses (11.8), (ll.9)ancf (II. 1 3)j (i = 2) hold. Then, the 

solution {u ,p) of problem ( 11 . 1 ) satisfies the following estimates 


(11.15) llulL. <a: 


+ (1 + a:’if) ||gi 


(11.16) IIpII. < (Ua;’if)||f|lx; + §i’§;'x(i^oc;'x)llg|lM: 


-1 


10-1. 


II, 2. The finite dimensional apDroxImation. 

We want to approximate problem (11.1) by a finite dimensional system. To this end, we 
introduce a discretization parameter 6 > 0 , and we assume we are given closed subspaces Xj^ and 
hjg (i = 1 , 2) contained in Xj and Mj respectively. The continuous bilinear forms a and bj are 
approximated by three continuous bilinear forms Og : X 2 g X Xi5 -> [R and bjj : Xjg X M ->■ IR 
(i= 1 , 2 ). 


For any f^ in X) and g^ in , we consider the following approximation of problem 
(11.1) -.Find (Uj , Pj) in X 25 X that 

(11.17) VVg€X^j, 9j('^S I ’ 

Vq5€M25, b2,(Uj .q^) = <g^ ,qg> . 

We assume that the forms a^ and bj^ satisfy the necessary and sufficient conditions for problem 
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(11.17) to be well-posed, as stated in Theorem 11.1 . More preciselg, for any g in (i = 1 , 2), we 
introduce the affine subspace 


(11.18) Kjg(g) ={v5€Xi5 ;Vqg€Mi, , b^^ (v^.qj = <g,q5> ) , 
and we set Kjj = Kjg(O). We make the following assumptions : 
a) there exists a constant > 0 such that 


(11.19) Vu,€K 


28 


sup 




^ (X 


'8"X 


18 ll'^8llX, 


(11.20) VVj€K, 5 \{ 0 } , sup ag(Ug,V 5 ) >0 ; 

^8 ^ *^28 

this second condition, if K,^and Kgg are finite dimensional, can be equivalently replaced by 

(11.21) dim K,^ = dim Kje, ; 

b) there exists a constant > 0 ( i = 1,2) such that 


(11.22), Vq,eM„. sup 
Moreover , we denote by 


^8^Xi8 


llvsllx.llq8lL 


^ &i8 > 0 . 


(11.23) ^ 5 = sup 

*^8 ^ ^28 ’ '^8 ^ ^18 


9r(Ur -Vr) 


■^6 "X, 


'8"X, 


the norm of the form a^ . 


Oornllaru 11.2 -.Assume that hupotheses (11.19), (11.20) and (11.22)j(i= 1,2) hold. Then 
problem (11.17) has a unique solution (Ug ,Pg). Moreover the solution (Ug ,Pg) satisfies the 
following estimates 


(11.24) II Ug lljj^ < *^18ll^8llx‘, ^28 ^ ^ *^18)(8^ II S^IIm^ , 

(11.25) IIPjIIm^ ^ ^18 ^ ^ ‘^18)(8^ ll^8llxi ^18 ^28 ^8 ^ ^ ^18^8^ II98 IIm2 


11.5. Error estimates. 

We assume now that there exist a solution (u,p) to problem (II. 1 ) and a solution (Ug , Pg) to 
problem (II. 1 7). We want to derive a general error estimate between them. 


First, we estimate how an element of Kj(g) can be approximated by an element of Ki 8 (Q 5 ). for 
any g and gg given in MJ . 

Proposition II. 1 ; Assume that hypothesis ( 1 1 .22)j holds. For any element v in Kj(g) , the 
following estimate is satisfied : 
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(11.26) 


inf llv-wjlx, 
W5 e K.s(g,) 


inf { IIV-V 5 IIX.+ sup 


(brM(Vs-qs) ^ 

n _ M J 


'^8 ^ ^18 


Q8 € M., 


8 I'M. 


sup 
^8 € 


< 9-98 ■ Q 8 > 

ihsIlM. 


] . 


Proof : Let Vg be any element in . By (11.22)) and [Br, Thm 0.1 ], there exists in such 
that 

(11.27) VqgeMjg. b),(z 5 . q^) = b) 5 (Vs , q^) - <05 . Q 5 > . 

(11.28) ||z,L< 6 u 5up II ’ ' 

' Ms«M| 

Clearly, the element Wg = Vg- Zg belongs to Kjg(gg) and satisfies 

(11.29) llv-Wgllx.^llv-Vgllx.dlZsISx) • 

Thanks to (11.27), we write 

b)g(Zg ,qg) = -bi(v-Vj ,qg) -(b)-b)g)(Vg ,pg) +<g-gg ,qg> ; 

hence, using (11.28) and the continuity of bj yields 


The last inequality, together with (11.29), gives the proposition. 


(b-bi8)(v8.98) , su, <g:9s - 98> ) 


+ sup 
98 € ^i8 


8 I'M; 


Next, we derive the main error estimate for u-Ug . 

Thenrem 11.2 : Assume that hypothesis (11.1 9) holds. Then the solutions (u,p) of (11.1 )and 
(Ug , Pj) of (11.17), satisfy the following estimate 


(11.30) 


llu-UgL ^c(Uai^)[(UXs) inj l|u-Wg|L 

2 Wge<2g(gg) 


+ inf 
Vg€X 


28 


{( 1 +K 8 )l|u-Vglix + sup 

^ 2g€X), 


(a- 3 g)(v 8 ,Zg) 



} 


+ inf {IIp-QsIIm + 
98^^18 


sup 

^8^ 


(b) b)g)(Zg , Pg) 

X 18 


) 


+ sup 

28^X15 


<f-^8 ■ ^8> 

11^8 ix, 


] 
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Proof : Let Wg be ang element in K 25 (gg). By (1 1. 1 9) , we have 

, X „ I. -1 ,z,) 

(II.31) IIU 5 -W 5 L sup ^ ^ ^ 


^8 ^ *^16 


-8"X 


But, thanks to ( 11 . 1 ), (11.17) and (11.18), for any in , we can write 

, Zg) + <f^ , Zj> 

= a(u, Zg) - . Z5) + b,(z3 . p) - <f-fg , Zg> 

Now, let Vg be any element in X 25 and be any element in . By (1 1 . 1 8 ) , we obtain 

- 25) = aCu.z^) - . z^) + b,(Zj , p-q^) + (bi-bi,)(z, . q,) - <f-fg , z,> 

hence. 


(11.32) 


as^Ug-Wg , 2 g) = a(u-Vg . z^) ^ a 3 (v^-W 3 , z,) + (a-a 3 )(v^ , z^) + b/z^ . p-q^) 


+ (b,-b,g)(z 3 .q^) - <f-f 5 ,z,> . 


Using (11.23) and the continuity of a and b, in the last formula, we deduce 
IIV'^sllxo ^ “u [ ^8 llu-wJL { (Y^-Xb) IIu-vJIx„ + sup 


(a-a,)(v, ,z,) 




" 8 llX, 


f ti 1 . ^b,-b,J(z^ , qj - <f-f. , Zf.> , 

+ {clIp-qsIlM, ^ sup L Y" ^ ^ sup ib^ ] - 


’ 2 ,€X ,3 


Zg€Xi3 I|2 bIIXi 


which gives the theorem. 


Remark 11.2 : Sometimes it can be useful to make a further assumption on the form b^g , which 
will always be satisfied in the applications we have in mind. It will allow us to prediot an optimal 
rate of decay for the error u-Ug even in the case where Mjg is a bad approximation of . Thus we 
assume that there exists a subspace 1H,gOf M,, such that the form b,g is in fact defined on 
X 28 X ^^18 satisfies 

(11.33) VVg€K,g,Vqg€H,g, b,g(Vg,qg)=0 . 

(In applications, H^g will be a subspace of M, containing M^g and such that the elements of M., can 

be approximated in an optimal way by elements of H,g). Then, we can replace inf by 

1^18 

inf in the right-hand side of (11.30). 

%€^18 

Now, we indicate a remarkable case in which estimate (11.30) can be simplified. The result 
follows easily from (11.32). 

Corollaru 11.5 : Assume that hypothesis (II. l 9) holds and that 
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(11.34) KijCK, . 

Then the solutions (u.p)of (ll.l)and (Ug,Pg)of {\\M) , satisfij the following estimate 


(11.35) 


inf Ilu-Wgilv 

2 Wj € K2j(gs) 2 


+ inf {(u-)f^)||u-Vg||x + sup 
'^8 ^ ^28 


(a-aj(v, , zj 




+ sup 
z^eX 


IS 


<f-h ■ 2 r> 
l|2sllx, 


] 


} 


Finally, an error estimate for p-pg is provided by the following result. 

Theorem 11.3 : Assume that huootheses (11.19)and {\\.22)^hold Then the solutions (u,p)of 
(II. Dancf (Ug , Pg)of (11.17), satisfy the following estimate 


Hp-PsHm inf , llu-wjlx, 

’ Wg € K2g(gg) 2 


(11,36) 


^ inf {0+Y8)llu-V8lL 
V8^^28 ' 


sup 

28«Xi8 


(a-ag)(Vg , Zg) 

ll^sllx, 


} 


+ inf {IIp-QsIIm + sup 
Ps^f^is 2g6X,g 


(bi-Pis)(2s.q8) . 

Il^sllx, ^ 


+ 


sup 

2 s^X,g 


<f-fs . ^8> 
lUsllx, 


] 


Proof : Let Pg be any element in M,g . By (11.22), , we have 

(11.37) < eu ■ 

’ Zg4X,g ll^sllx, 

But, thanks to problems (11.1) and (11.17), for any Zg inX,gand Vg in X 2 g we can write 

1 Ps~Q 8) - ~ Sg(Ug , Zg) + <fg , Zg> — hijCZg , Pg) 

= a(u, Zg) - 8g(Ug , Zg) + b,(Zg . p) - <f-fg . Zg> - b,g(Zg . Pg) 

= a(u-Vg , Zg) - 3g(Ug-Vg , Zg) + (a-ag)(Vg , Zg) + b,(Zg , p-Pg) 

^-(b,-b,g)(Zg,Pg)-<f-fg,Z5> 

Using (11.23) and the continuity of a and b, in the last formula, we deduce from (11.37) 
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b-^sIIm ^5 llu-Usilx +{(^+Xs)liu-vJlx + sup ’ V 




^5 NX, 


f II II , qj , <f-f,. , z,.> 

+ {cIIP-QsIIm + sup ^ ^} + sup ^^-S- 




^NX, 


2c € X 


IS 


H NX, 


This, together with (11.30), proves the theorem. 


Remark 11.5 : Assume that hypotheses (11.19) and ( 1 1 .22) j hold and that the discrete forms a^ 
and bjj (i = 1,2) are the restrictions to the spaces of discretization of the continuous forms a and 
bj (i = l, 2) respectively, and that f and f^ coincide. Then, the previous estimates can be 
substancially simplified as follows : 

(11.38) ||u-uJL < c(l+a'^) inf lIu-w^lL , 

' w, € K2s(g,) 2 

and 

(11.39) IIp-p^IIm < c(l+g;s)(Ua'j) ( inf Hu-wjlx + inf IIp-QsIIm, ) ■ 

' w,€K2s(gs) 2 ’ 


III. Variational formulation of the Stokes problem in weighted 
Snhniftv spaces. 


We are interested in solving the Stokes problem in the domain Q = ]-1 ,1 [ : thus, given a 
force field f in Q and a viscosity v > 0, we look for a velocity field u and a pressure p (defined up 
to an additive constant) such that the following equations are satisfied 


(lll.l) 


-vAu+gradp=f in 0 


[ divu = 0 in Q , 

together with a Dirichlet boundary condition on the boundary 90 of 0; we shall first study the 

homogeneous case 

(III. 2) u = 0 on 90 

Since we want to study a spectral Chebyshev approximation of the Stokes problem , we first give a 
variational formulation of (III.1)(I11.2) in terms of weighted Sobolev spaces, the weight being 
precisely the Chebyshev one in each variable. Then we use Theorem 1 1.1 to prove the well- 
posedness of the variational problem. Finally, we extend the results to the case of non homogeneous 
boundary conditions. 


111.1. The weighted spaces and the homogeneous Stokes problem. 

Let us briefly recall some basic material about weighted spaces of Chebyshev type (for 
further details, see e.g. [CHQZ][Ma2]). If q( 0 = ( 1 denotes the Chebyshev weight on the 

interval ]-l ,1 [, let 

Lj(-l .1 ) = { >P : ]-M [ ^ IR ; lo D^(0 e(C) dC < -Hoo } 
be the Lebesgue space associated with the measure q(C) dt; , provided with the inner product 
(III. 3) (<P.'l»)g = lo ^(C) 'I'(C) Q(C) dC 

and the norm IMIq^ = (.,.)J^^ 

A scale of weighted Sobolev spaces is defined as follows : for any integer m > 0, H^(-1 ,1) 
is the subspace of L^(-l ,1) of the functions such that their distributional derivatives of order 
< m belong to Lp(-1 ,0 ; it is a Hilbert space for the inner product associated with the norm 
(UI.4) = . 

where 

(III. 5) |(plk„=||d^/dc'|lo,p • 

For a real number s = m + a,0<a< 1, we define H®(-1 ,1) as the interpolation space between 
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H^(-l ,1) and ,1 ) of index a (cf. [LM]); we denote its norm by 1|.1|^^ . 

Finally, we can apply a rotation and a translation to define similar Sobolev spaces on any 
segment of length 2 in We use the same notations as before to indicate them , as well as their 
norms. 

The generic point in the square Q will be denoted by x = (x.y). We introduce the vertices 
8j , J € Z/4Z, of Q (where aj^, follows Sj counterclockwise) , and call Fj the edge with vertices 
3j_, andaj ; for any edge Fj , J € Z/4Z, iij is the unit outward normal to Q on Fj and ij the unit 
vector orthogonal to nj , directed counterclockwise. 



Figure III. 1 
The square Q. 


The Chebyshev weight on 0 is defined as co(x) = g(x) p(y). Let 
L^(0) = { V : Q -»• [R ; ifi^(x) w(x) dx < +oo } 
be the Lebesgue space associated with the measure gj(x) dx, provided with the inner product 

(111.6) ('P>'F)oj = 1o >P(x) >V(x) w(x) dx 
and the norm IMIo,^ = (.. )^^^ 

Next, a scale of weighted Sobolev spaces is defined as follows : for any integer m > 0, 
H|^(Q) is the subspace of L^(Q) of the functions such that their distributional derivatives of 
order ^ m belong to L^_^(0) ; it is a Hilbert space for the inner product associated with the norm 

(111.7) lklU,„ = (Lr=o . 

where 
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For a real number s = m+oc,0<a<l,we define H^(Q) as the interpolation space between 
H'^(Q) and of index a; we denote its norm by ||.|| 

Being concerned with homogeneous Dirichlet boundary conditions, for any integer m ^ 1 . 
we consider the closed subspace of the functions of H|^(Q) which vanish on the boundary 30 
together with all their derivatives of order up to m-1 (the traces being defined in the sense of 
[LM]) ; this space, denoted by H[^q(Q), is the closure of DiCl) under the norm of H^(Q) (see 
[BM, Prop. 11.9]). Due to the Poincare inequality, an equivalent norm on H|^g(0) is the seminorm 
. The dual space of H|])q(Q) will be denoted by H'""(Q) ; whenever the space L^(Q) is 
identified to its dual space, we have for instance 

(111. 9) H-^(Q) = {f + 3g/8x + 8h/8y,(f,g,h) € [L^(Q)]^} . 

Now, we go back to the Stokes problem (lll.1)(lll.2). Assume that f belongs to [H"\q)]^. A 
weighted variational formulation of (III.1) is obtained by requiring that the first equation in 

(1 1 1.1) is satisfied in [H’’(Q)]^, and the second equation in L^(Q). In order to satisfy also 

(1 1 1.2) , we define the space 

(111.10) X = X, =X2 = [H^g(Q)]2 , 

in which we look for the velocity. Since the pressure p is defined up to an additive constant and 
div u has zero average in Q , we introduce the closed subspaces 

(111.11) M, = {qe L^(Q) ; q(x) cj(x) dx = 0 } 
and 

(111.12) Mg = { q € L^(Q) ; q(x) dx = 0 } . 

in which respectively we look for the pressure and we choose the test functions to enforce the 
divergence-free condition ; this choice will be justified later. 

Thus, for any f in X', we consider the following variational formulation of the Stokes 
problem (111. 1 )(l 11.2) ; f/nof (u,p) X x ^ such that 

(111.13) VveX, V grad u . grad (vco) dx - div (vw) p dx = <f,v> , 

V q € Mg , div u q cj dx = 0 

This formulation will turn out to be a particular case of the abstract variational system (II. 1). 
provided we define the bilinear forms a : X x X IR and bj : X x M, ^ [R (i = 1,2) respectively 
by 

(111.14) a(u,v) =- V <Au,v> = V grad u . grad (vco) dx , 

(111.15) b,(v,q) = <v. grad q > = - div (vto) q dx , 
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(111.16) b 2 (v,q) = - (d1vv,q)^^ = - divvq CO dx . 

Note that we have 

(111.17) K, = {v€ X ;div (vco) = 0 in Q } 
and 

(111.18) <2 = { V 6 X ; div V = 0 in 0 } . 

We have at once the following result. 

Proposition III. 1 : The forms a : X x X ^ [R and bj : X x Mj -> DR (i = 1 , 2) are continuous. 


Proof : It is a straightforward consequence of the continuity of the mapping 
ip -> co"^ grad ((pco) from H^q(Q) into L^(0). In order to prove this result, we write 


(0 ’ grad (ipco) = grad (p + ip 


f x/(1-x2)) 
y/(1-y^) 


and, due to Hardy's inequality (see [N, Chap. 6, Lemme 2.1]), the terms ip x/(l-x^) and 
ip y/( 1 -y^) can be bounded by the norm of ip in H^q(O) (see [CQ2, Lemma 1 ,2] for details). 


Thus, in order to apply Theorem 11.1 to problem (111.13), we have to check the inf-sup 
conditions (ll.8)(ll,9) and (II, 1 3)j (i = 1,2). 


III. 2. The inf-suD condition for a. 

Let us first deal with the form a. For some technical reasons, we introduce the weighted 
Sobolev spaces relative to the inverse of the Chebyshev weight : for any real number s > 0, the 
spaces H®/p(-1 ,1 ) are defined in the same way as the spaces H®(-1 ,1 ) with p replaced by 1 /p and 
provided with the norm ll llsi/p ; the spaces H®/^^(Q) are defined in the same way as the spaces 
H^(0) with CO replaced by 1/to and provided with the norm ||.||g ; for any integer m > 0, we 

denote by H'J’/j^q( 0) the closure of £>(0) under the norm 11.11^^ . We recall the following result 

due to [BM], 

Lemma 111.1 ; For any integer m > 0, the mapping : (p -*■ to^^^ ip is an isomorphism from 
H[|^ o(Q) onto H^(Q) and from H^(Q) onto H'j"/^o(Q) . 

Now, we are in a position to prove the following proposition. 

Proposition III. 2 : There exists a constant a > 0 such that 

(111,19) Vu€Kj, sup -5^ s a Hull, „ , 

V 6 K, 
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(111.20) Vv€K, , sup ^a||v||,^ . 

u € Kg 

Proof : Let us start with (III. 1 9). If u belongs to Kg , then there exists (p in Hq(Q) such that 
u = curl (p (see e.g. [OR, Chapter 1 , Thm 3. 1 ]); since u is in X, ip belongs in fact to H^q(Q). 
Definev = co"' curl (ipu). By Lemma lll.l , (pw belongs to H^/^^q(Q), hencev belongs to H^q(Q), 
with 

l|v||, ^^ <cl|u||,^ • 

Moreover, 

a(u,v) = V (curl u) (curl (vcj)) dx = v Aip A(ipw) dx 
According to [MM , Lemma 3.2] , there exists a constant c> 0 such that 
Jo A4) A(ipco) dx >c||(p|lg^ . 

We conclude that (111.19) holds. 

In order to check (111.20), we use a similar argument. If v belongs to Kg , then there exists 
ip in q(Q) such that vco = curl ip. We set u = curl (ipco”’), so that u belongs to K, , and we 
conclude as before. 

111.3, The inf-suD condition for b - (i = 1. 2). 

Next, we check the inf-sup condition (11. 1 3)j for the forms bj (i = 1 , 2). To this end, we 
recall some trace and regularity results in weighted Sobolev spaces, due to [BMj. 

With the notations of Figure 111. 1 , we have the following lemmas. 

Lemma III. 2 [BM,Thm 11.3]: The trace operator R : v -> ( vip )j ^2/42 is linear continuous 

from the space hJJ(Q) onto the subspace of IT j ez/42 functions ( »Pj )j ^ 2/42 

satisfying 

(111.21) '^J^7L/47L, (pj(aj) = (pj^,(aj) . 

Moreover, the operator R admits a continuous right inverse. 

Lemma 111.3 [BM, Thm 11.4] : The trace operator S : v -*■ { V|f^ , 9v/9nj )j ^2/42 tinear 
continuous from the space H^(Q) onto the subspace of rTjg2/42 ^ 

functions ( tpj . Xj )j€2/42 satisfying 
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(111.22) VJ6Z/4Z, I ipj(aj) = 

(a(py8lj)(aj) = Xj,,(aj) . 

Moreover, the operators admits a continuous right inverse. 


Lemma III. 4 [BM, Thm 11.2] : The trace operator T : v -> ( V|p^ , av/arij )j ^ 2/42 Unear 

continuous from the space onto the subspace of TTjg2/42 x of 

the functions ipj . )j ^ OH. 21). Moreover, the operator T admits a 

continuous right inverse. 


Finallg, consider the Dirichlet problem 


(111.23) 


- Alp = X in Q 
Ip = 0 on ao 


Lemma III. 5 [BM, Thm 111.1] -.If x belongs to LmQ), then the solution ip belongs to 


H^(Q) riH’(Q) and 


(111.24) lkll2..<c||x||o^^ . 

If X belongs to L^/^^(Q), then the solution ip belongs to H^/^(Q) Pl Hq(0) and 

(111.25) ll'flll2,1/to ^ Iloilo, l/cj 


Now, we are in a position to prove the following propositions. 
Proposition III.5 ; There exists a constant ^2 > 0 

(111.26) VqeM 2, sup ^& 2llqllou ■ 

v€X llvlli.u, 


Proof : Let q be any function in M2 . We consider the unique solution ip of (111.23) with x = Q. for 
which we have 

Iloilo, u • 

We want to find a function \p in H^^(Q) such that 

(111.27) curl >p = - grad ip on ao 
or equivalently 

(111.28) '^S^TU^TL, aip/aij = - aip/anj and dip/anj = dip/dij = 0 on fj 
To this end, we set 

(111.29) iP|(a,)=0 , 
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and, for J = II . Ill , IV, I successively. 

(111.30) Vx€Tj, H(j(x) = - fpj (9(p/anj) x([aj_, ,x]) dcr , 

where x([aj_, ,x]) is the characteristic function of [aj_, ,x] in Tj . By Lemma 111.3, 9(p/9nj 
belongs to and satisfies 

(111.31) VJeZ/4Z, [ (a(p/3nj^,)(aj) =(3(])/aij)(aj) =0 , 

(a<p/anj)(aj) = - (a<p/aTj_^,)(aj) = 0 . 

Hence, ii/j is in Lg(fj) and aij/j/aij = - aip/aoj is in H^^'’(rj), so that belongs to hJ^'’(Tj) 
(see [BM. Lemma 11.4]). Moreover, since q is in Mg , we have 

^je 2 Z /42 Ifj (9^/9rij) dcr = A(pdx=-|jj qdx = 0 , 

hence, the functions H(j satisfy 

VJ€Z/4Z, H'jCaj) = <t<j^,(aj) 

We also have by (111.31 ) 

VJ€Z/4Z, (ai|>yaTj)(aj.,) =(aHfj/8Tj)(aj) = 0 . 

Thus, the functions ( , 0 )jg 2/42 satisfy the compatibility conditions (111,22). By Lemma 

111.3, there exists a function vjj in H^(0) such that 

»!; = H/j and a»|;/anj = 0 on fj . J € 2/4Z , 

and which satisfies 

(111.32) Il'l*ll2.w ^ ^J€Z/4Z ^ 

Finally, we set v = grad ip + curl i|; . Thanks to (111.27) and (111.32), the function v belongs to X 
and satisfies 

(111.33) llvl|,^<cllq|lo^ . 

Moreover , we have div v = q, so that 

b 2 (v,q) = q^ cj(x) dx , 

This formula, together with (111.33), gives the proposition. 

Remark III. 1 : Note that the proof of Proposition 1 1 1.2 utilizes the condition q(x) dx = 0 in 
the definition (111.12) of Mg . 

Proposition III. 4 : There exists a constant g, > 0 such that 

(111.34) VqeM, , sup Halloo. • 

v€X ll''lli.w 

Proof : The proof is similar to the previous one. Given any function q in M, , denote by ip the 
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solution of the Dirichlet problem (111.23) with y = qw. Since qu is in Lemma Hi. 5 

implies that tp belongs to H?/jQ).with 
Il'Plki/u^cIlqllo^ ■ 

Next, we build up a function in such that again (111.27) holds. First, we define its 

trace on 90 by (111.29) and (111.30). By Lemma 111.4, 8ip/9nj belongs to Hj^^(fj) e L’(rj), 
thus »j/j is well-defined and belongs to H^/p(rj). Moreover, since q is in M, , ipj satisfies 
VJ€Z/4Z, 'Pj(aj) = »1'j+i(aj) 

Applying again Lemma 1 1 1.4, we see that there exists a function ip in H^^^_^(Q) satisfying 
ip = ip j and Sip/anj = 0 on fj , J € Z/4Z 
and such that 

(111.35) II*FII2,1/u> ^ £ 2/42 Il'^’jll5/4,1/p ^ ^ 

The function v = grad tp + curl tp belongs to X and satisfies 

(111.36) llvll, ,/^<c||qllo^ , 

Finally, we set V = vco"V By Lemma 111.1, 

Wi.u, <c||v||, ■ 

while on the other hand div (vco) = div v = - qco , so that 
b,(v, q) - [fj q^ cj(x) dx . 

This proves the proposition. 

Remark 111.2 ; Here, we have used the condition q(x) cj(x) dx = 0 in the definition (111.1 1 ) 
of M, . 

111.4. The existence and uniqueness results. 

Thanks to Propositions III. 2, III. 3 and III. 4, we can apply Theorem 11.1 to the variational 
problem (111.13), and obtain the main result of this section. 

Theorem III.1 : For each f in X' , there exists a unique variational solution (u,p) in X x M, to 
the Stokes problem (III. 1 )(lll.2). Moreover, the following inequality is satisfied for a constant 
c > 0 

(111.37) llul|,^-H|lp||o^<c||f||x. . 

Remark III. 3 : In the Stokes problem (111. 1 )(lll.2), the pressure is defined up to an additive 
constant. In the variational problem (111.13), this constant is fixed by the condition 

(111.38) jfj p(x) (j(x) dx = 0 
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(which is not the usual one - see [GR, Chapter I, Thm 5.1 ][BMM2, (V.2)] for instance). 

We conclude this section with the case of non homogeneous Dirichlet boundary conditions, 
which we write in the form 

(111.39) u = tpj ouFj, 0€l/47L . 

The following result holds. 

Theorem III.2 : For each f in X' and for each ((pj)jj 2/42 in njez/ 4 z satisfying 

(111.40) yd^TL/4TL, ipj(aj) = ipj^i(aj) 

(111.41) Z:j,z/4Z Jrj = 0 ’ 

there exists a unique variational solution (u,p) in [H^(Q)]^ x M , fo the Stokes problem 
(lll.l)(lll.39). Moreover, the following inequality is satisfied for a constant c> 0 

(111.42) l|ulli,u> IIpIIo.u) ^ “^ ( 11^ llx' € 2/4Z Il'Pj Il3/4,e ) 

Proof ; In order to apply Theorem 111. 1 , we look for a function in H^(0) such that 

(111.43) divu[, = 0 in Q and U|, = ipj on Tj , J € Z/4Z 

or equivalently, setting u,, = curl i|;, for a function in H^(Q) such that 

(111.44) VJ€zZ/4Z, 94</5Ij = ipj . Hj and 0ii;/anj = - cpj .ij on Tj 
To this end, we set 

(111.45) >i«|(a|)=0 , 

and, for J = II , 111 , IV, I successively, 

(111.46) VxeTj, 4 »j(x) = »i<j.i(aj_,) + (ipj . rij) x([aj_, , x]) do , 

where x([aj_i , x]) is the characteristic function of [aj_, , x] in fj . For J € TL/ATL, the pair 
( 4 fj , - (pj .ij) belongs to Hj'^'*(rj) x H^^'^‘^(Fj) and, due to (111.40) and (111.41), we have the 
conditions 

VJ€Z/4Z, H'j(aj) = >pj*i(aj) , 

(3ipj/aTj)(aj) = -(ipj,i . nj^,) (aj) , 

(iPj .nj) (aj) = (ai;j^i/atj^i)(aj) . 

It follows from Lemma 111.3 that there exists in H^(Q) satisfying 
VJeZ/4Z, ip = ipj and aip/anj = - ipj . Tj on Fj , 
whence (111.44); moreover, we have 
lkll2,u, ^ LjsZ/ 42 Il4ljll3/4,p 

Next, we define Uj, = curl \p and, using Theorem lll.l , we consider the unique solution (u,p) of 
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problem (111.13) with <f,v> replaced by <f,v> - a(uj, ,v). Clearly, the pair (u = u+u^^ , p) is 
solution of the Stokes problem (1 1 1. 1 )(l 1 1.39) and satisfies (1 1 1.42). 

Remark III. 4 : The previous result can also be proven by a direct application of the abstract 
Theorem II. 1 . As a matter of fact, if ( iPj )j g 2/42 belongs to f[j ^2/42 and satisfies 

(111.40), by Lemma III. 2, there exists a vector field in [H^^(0)]^ such that 

(111.47) u;|r^=:4'j ,J€Z/4Z 

and 

(111.48) < C ZIjg2/42 Il'f’jll3/4,p 

Let (u' ,p) be the solution in X x of the variational problem 

(111.49) VveX, a(u' ,v) + b,(v,p) = <f,v> - a(u^ , v) , 

VqeM 2, b2(u’,q) = - b2(u^, q) . 

Thanks to Theorem II. 1 , the continuity of the forms a and bg and (111.48), we have 

l|u + IIpIIo.u ^ ^ ^ ^J€2/42 ll‘Pjll3/4,p ^ 

Finally, let us set u = u' + . Then, (u,p) is the solution of the Stokes problem (II 1. 1 )( 1 1 1.39). 

Indeed, the only nontrivial condition to check is the continuity equation. By (1 1 1.49) we have 

(111.50) V q€ L^(0) / q(x) dx = 0, div u q co dx = 0 . 

Condition (111.41) implies that 

Iq div u dx = 0 , 

which, together with (111.50), yields div u = 0 in Q. 
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IV. A Galerkin method for the Stokes problem. 

Henceforward, we fix an integer N ^ 3. In this section, we will study the convergence 
properties of a spectral Galerkin approximation to problem (111.13). Such a method is never used 
in practical computations - being less efficient than a tau or a collocation scheme. However, we 
consider it here because it is the simplest method which can be handled by the theory of Section II 
and in which the main difficulties are present. 

IV. 1. The discrete problem. 

We look for an approximate solution of (111.13), the components of which are algebraic 
polynomials of degree ^ N in each variable. This solution will be defined by a Galerkin method 
using the variational formulation of Section III. A suitable orthogonal basis in order to study the 
algorithm consists of the Chebyshev polynomials of the first kind. 

Let us begin by recalling some basic material about Chebyshev expansions. We denote by 
T^(x), m = 0, 1 , ... the Chebyshev polynomials of the first kind. They are defined by 
T^(x) = cos (m Arcos x) , and satisfy the orthogonality relation 
CV'l) TjC:)T„(C)(l-C^)‘’'^dC = c^(TT/2)6^„ ,m>0,n>0 , 

where Cq is equal to 2 and c^ is equal to 1 for m > 1 , and 6^^^, denotes the Kronecker symbol. 
Moreover, for any m ^ 1 , the Chebyshev polynomials satisfy the differential equation 
(IV.2) VC€]-1,1[, ((l-C^)’^^OC))' + m^TjC)(1-C;^)'’^^ = 0 

and the relations 

(IV.3) VC€]-1.1[, T^.,(C) = 2CTJ0-T,.i(0 . 

(IV.4) VC€]-1,1[, TjC) = T;„,,(i;)/2(m+1)-T;„.,(C;)/2(m-1) . 

Let Pfj(- 1 , 1 ) denote the space of the algebraic polynomials of degree < N in one variable, 
restricted to the interval ]-1,1[. Each ip in Pn(- 1,1) can be expanded as 
ipCC) = L1 o ^PmTJO.with 

(IV.5) ip^ = (2/rccj|', 4«(C)UC)(1-C^)'’'^dC ,0<m<N . 

Pn(- 1 ,1 ) will be the subspace Pn(-1 .OH Hq(-1 ,1 ) of the polynomials vanishing at the end 
points C; = ± 1 • 

Next, we denote by .0]®^ the space of the algebraic polynomials in IR^ 

which are of degree < N in each variable. Finally, we set PN(f^) = Hq(Q). 
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Let us now introduce the Galerkin approximation to problem (111.13). We look for the 
approximate velocity in the space = [Pn(^))]^ and for the approximate pressure p*^ in 
P^(Q). Weconsider the following problem : Find (u^, p^) in x P,^(0) such that 
(IV. 6) Vv€X^, V grad . grad (vco) dx - jfj div (vw) p'^ dx = <f,v> . 

Vq€P|^(0), div q cj dx = 0 . 

However, as it will be discussed in the next subsection, in order to have a well-posed discrete 
problem, we must restrict the space of the pressures to a proper subspace of P^(Q). 
Similarly we restrict the space of test functions for the divergence free-condition to a proper 
subspace M 2 N of P^(0). Then, we obtain a particular case of the abstract approximate problem 
(II. 1 7), if we set X ,5 = Xgg = X|,j provided with the norm of X, and Mgg = Mgjj provided 

with the norm of L^(Q), and if the forms a^ and b^^ (i = 1 , 2) are respectively the forms a and h- 
(i = 1 , 2) defined in (111.14) to (111.16). 

In order to apply the abstract convergence results of Section 1 1 , we have to choose the finite 
dimensional spaces and for the pressure, in such a way that the inf-sup conditions for 
the above forms hold. In the next subsection, we characterize those pressures which cannot satisfy 
such a condition for the forms b, and bg . 


IV. 2. The spurious modes of the pressure. 

Spurious or "parasitic" modes of the pressure are those components of the numerical 
pressure that are not controlled by the discrete equations which approximate the Stokes system. 
Parasitic modes in finite difference or finite element methods have longly been investigated (see 
e.g. [GR]). In spectral methods, attention on such a problem was first brought by Y. MORCHOISNE 
[Mo]. The characterization of spurious modes for various spectral methods of mixed 
Fourier-Legendreor fully Legendre type has been carried out by [BMM 1 ][BMM2]. 


Our aim is to characterize the polynomials q in P^(f^) which satisfy, respectively for i = 1 
or 2, the following condition 
(IV.7)j VveX^, bj(v,q)=0 . 

Let Zjfj denote the subspace of all q in Pjj(O) for which (IV.7)j holds. It is clear that such 
polynomials cannot verify (II. 1 3); . Moreover, if M,^ contains a non trivial element of , the 
pressure in the solution of the discrete problem is not unique. 
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Let us first deal with the form b, defined by (III. 1 5). We recall that the polynomials 

m = 1 N+ 1 , form a basis of P^C- 1 .1) which, due to (IV. 1 ) and (I V.2) , is orthogonal with 

respect to the inner product 

(IV.8) (u.v)i/^ = I^, u(c;)v(c;)(l-C;^)'^^dC . 

Precisely, one has 

(IV.9) l’, T;j:c) T;,(0 (l-C^)'^^dC; = (Ti/2) m^8^^ ,m^0,n>0 . 

Lemma IV. 1 : The dimension of the subspace is equal to 8. 

N 

Proof ; Let us expand each q in as q(x) = L^^nzO ^mn Condition (IV. 7), is 

equivalent to the conditions 

(IV. 10) Vv6P°(Q), (0(vco)/ax) qdx = 0 

and 

(IV. 11) Vv€P;(0), jo (3(vco)/ay)qdx = 0 . 

Let us consider first (IV. 1 0). A basis in Pn( 0) is given by the polynomials 

(1V.12) v^„(x) = (l-x^)T;„(x)[T„(y)-T„(^,(y)] , 1 < m ^ N- 1 , 2 < n ^ N , 

where oc(n) is equal to 0 if n is even and to 1 if n is odd (we use here the fact that 

T^(± 1 ) = (± 1 )"’). From equation (IV. 2) and the orthogonality condition (IV. 1 ) , it is clear that 

(IV, 1 0) is equivalent to the set of relations 

CV.13) q^n - c„(„) = 0 , 1 ^ m < N-1 , 2 ^ n < N . 

Working out condition (IV. 1 1) in a perfectly symmetric way, we end up with another set of 
relations, namely 

CV H) q^n - y«(m)n = 0 . 2 < m ^ N , 1 < H ^ N- 1 . 

The relations (IV. 1 3) and (IV. 1 4) provide an orthogonal basis for 2,,,^ . This is given by the 
four modes TQ(x)Tg(y), T,,j(x)TQ(y) , TQ(x)T,g(y) and T,,j(x)T|^(y) (since the corresponding 
coefficients do not appear in (IV. 1 3) and (IV. 1 4)), and by four other polynomials, the non-zero 
Chebyshev coefficients of which are respectively Po, , q,Q , q, , and Pqj plus the coefficients defined 
from these by the relations (IV. 1 3) and (IV. 1 4). This proves the lemma. 

The next lemma allows us to find another basis for , which is easier to handle. 

Lemme IV. 2 ; The set of all elements in P^^- ' ) which are orthogonal with respect to the 

inner product i.,.)^to the space Pn(-1 .1) is the subspace of dimension 2 spanned by 

^ ■ ^N+1 } ■ 
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Proof : Each polynomial ip in P,g(-l,l) satisfying ip(±1) = 0 can be written as 
ip(0 = (1 -C^) ip(C) , where ip belongs to Pn_ 2 (- 1 .D- Expanding ip according to the basis 
( <m<N-i (IV. 9) , we see that ip is orthogonal to and . On the other hand, the 

subspace P|^(-l,1) = { tp e Pf^(-1,l) ; ip(±1) = 0 } has codimension 2, hence the lemma is 
proved. 

Cnrnllaru IV. 1 : The set of all elements in P^(0) which are orthogonal with respect to the 
inner product i■,■)^^to the space { v 6 P|g(0)) ; v(± 1 ,± 1 ) = 0 } /s the subspace of dimension 4 
spanned by . 

Proof : Each polynomial v in P^COl) satisfying v(+ 1 ,± 1 ) = 0 can be written as v = w + z, with 
V (x.y) € n, w(±1,y) = z(x,±l) =0 

(take for instance w(x,y) = (( 1 +y)/2) v(x,l) + ((1 -y)/2) v(x,-l ) and z = v - w). Hence, it 
follows from Lemma IV. 2 that ( Tilj , is orthogonal to { v € P,g(Q) ; v(+ 1 ,± 1 ) = 0 }. On 

the other hand, since this space has codimension 4, the corollary is proved. 

Noting that, if v belongs to X|,j , then w"’ div (vio) belongs to P^(0)) and is equal to 0 in 
(± 1 .+ 1 ) , and recalling the proof of Lemma IV. 1 , we obtain the following characterization. 
Proposition IV. 1 : The subset is the vector space of dimension 8 spanned by { Tq , 
and by )®^ 

Remark IV. 1 : Note that the basis in { T,;, , )®^ is a simple linear combination of the 

orthogonal basis of Z,,g defined in the proof of Lemma IV. 1 , as it can be seen by using (IV. 4). 

We consider now the form b 2 defined by (III. 1 0). First, we recall the following result, due 
to [BMM2, Corollaire V.1]. Let us for a while denote by L^^ , m = 0, 1, ... the Legendre 
polynomials. 

Lemma IV. 5 : The range of by the divergence operator is the subspace 0,,^ of P|,,(0) of 
codimension 8 defined by 

(IV. 15) D^ = { r € P^(0) ; r(±l .±1) rOand V qe { Lq . LJ®^ r(x) q(x) dx = 0 } , 


Remark IV. 2 ; The lemma can also be established by a proof similar to the one of Lemma IV. 2, 
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by using the properties of the Legendre polynomials. 

We can now characterize the space . Let us introduce the polynomial in P,^(-1 ,1 ) 
defined by 

(IV.16) V ip € P^(-l,l), 1^1 qN(C)<p(C)Q(C)dt; = f^ ipCOdt; 

(i.e. q^j is the orthogonal projection of 1/p onto P^^C-1 ,1) with respect to the inner product 

Proposition IV.2 : The subset ts the vector space of dimension 8 spanned by { q^^ , 
and by {T' ,T'^i }®^ 

Proof : The subspace is precisely the orthogonal space to D,^ in P^CO) with respect to the 
inner product (.,.)^ . It follows that has dimension 8. Next, from (IV, 1 5) and Lemma IV.2. we 
deduce that contains the subspace spanned by { Tj;, , }®^. In order to obtain the remaining 

components, we have to translate the orthogonality relation in (IV.16) in terms of the inner 
product (-,-),_j. To this purpose, we observe that, for any ip in Pfj(-1 ,1 ), if X is the coeficient of 
in r, we have (see [DR, § 1.13]) 

Li ip(UL„(i;)dl;=xL, =>'2“(NI)^/(2N)!(N»I/2) 

(♦,T„)„ = X(t;",V„ = Xn2-'' , 

so that 

(IV.17) V4)€Pn(- 1,D, j’, (p(OL^(t;)dt;=[2^(N!)^/(2N)!(N+l/2)]((p,T^)^ . 

Using this property and the definition (IV. 1 6) of , we deduce from (IV. 15) that Z 2 ^ contains the 
subspace spanned by { q^j , T^ }®^. 

Finally, let us check that the polynomials T,;, , T,;,^i , q^ and T^j are linearly independent 
(which will imply that the 8 elements of ( q,g , T^ )®^ U { T^ , T,;,^, }®^ are linearly independent). 
Assume that 

X, q|g + X 2 T^ + T^ + Pg = 0 

Using (IV. 4), we have 

Xi qN = (X2/2(N-1))T‘_, -Pi T' -(P2 +X2/2(N+1))T^^, . 

Since N is ^ 3. T’., = l is orthogonal to TJg_, , Tj|, and with respect to (.,.),/g , so that 
0 = X, j’, q^(0(l-t;^)’^^dC;=X, 1^, (1-t;^)dC =4X,/3 , 

whence X, = 0. Since the T,[ 's are mutually orthogonal with respect with (., )i/p , we obtain at 
once Xg = y, = P 2 = 0- 
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Remark IV. 5 : In this work, we are not concerned with the Stokes problem in the cube ]-l , 1 [^■ 
However, note that, in the three-dimensional case, the dimension of the space Z||,j (i = 1 , 2) 
depends on N : it is equal to 1 2N+4. Indeed, by the same techniques as in [BMM2, §5], it is an easy 
matter to prove that (resp. Zg^^) is spanned by { Tq , (resp. { q,,j , }®^) and by the 

( 1 2N-4) independent modes of 

U .1)®{ )©{ }] 


IV. 5. An inf-suD condition for th e forms b. (i = 1, 2^. 

The characterization carried out in the previous section suggests the most direct choice of 
the test and trial spaces and for the pressure, Precisely, we choose for M||.j the orthogonal 
complement to Zj,,j in P|g(Q) (i = 1 , 2) , i.e, 

(lV.18)j = { q e P^(Q) ; V r e Zjfj , (q,r)^ = 0) . 

Due to (IV.7)j , Mj,g is characterized as 
(IV. 19) = { CJ'^ div (vu), V € Xfj ) 

and 

(IV. 20) M2^ = {divv,v€X^} . 

Note also that is contained in Mj (i = 1 , 2). 

In order to check the inf-sup conditions (ll.22)j for bj over the spaces X,,j x (i = 1 , 2), 
we apply an abstract result of tensor algebra due to [Be, Chapter V, Appendix B], which we recall 
here. Let H be a Hilbert space, and denote by (,,.) its inner product and by ||.|| the corresponding 
norm. For any pair of planes A and B in H , we define the gap between A and B as the quantity 
(IV. 21) 6(A,B) = inf { ||ip-H;|| ; ip € A, 4; € B and||(|i|l = |1 h;1I = 1 } . 

Setting 

(IV. 22) i](A,B) = sup { (ip,»|;) ; (p 6 A, ip € B and l|(p|l = Ikll = 1 } , 

one has 

(IV. 23) 6^(A,B) = 2(1-t](A,B)) , 

Lemma IV. 4 : Let (A,B) be a pair of planes in H, the intersection of which is {0}. Any 
element q in H®^ which belongs to (A®^)'*’ and to (B®^)'*' can be written 
(IV, 24) q=r+s , 

where r and s belong to A‘*'©B‘*' and B'*‘©A'*' respectively and satisfy 
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(IV.25) sup {||rl| ,||s||} <(U2/8(A.B))^||q|| . 

Proposition IV.5 : Let the space be defined by (IV. 1 8)^ . There exists a constant |, > 0 
independent of N such that 

(1V.26) Vq€M,,. sup > li INlou • 

Proof : In the space H = P^(-1 .1) provided with the inner product t,.)^ , define A = { Tg , } 

and B = { T,;, , }, and recall that, due to Lemma IV. 2, B-^ is the space ( 1 -C^) Pn_2(-1 .0. 

According to Proposition IV. 1 , is equal to (A®^)‘‘‘ fl (B®^)'^. Since N is ^ 3, A fl B is {0}, 
hence, by the previous lemma, each q in can be written as in (IV. 24), with 
(IV.27) r(x,y) =(l-y^) LI’J a^(y)Tjx) . 

and 

(IV.28) s(x,y) = (1-x^) LI'J b„(x) Vy) . ^ € P^_ 2 (-l ,1) . 

Using (IV. 2), we write 

rco = ( 1 lJIjJ ajy) T Jx) ( 1 

= (8/ax)[( 1 ( 1 lI'I ajy) T;,(x)/m^] = a(vco)/8x , 

where V, defined by 

v(x) = - ( 1 -x^) ( 1 -y^) lI'J a^(y) i;„(x)/m^ , 

belongs to P,^(Q). Moreover, we have 

Ilav/8xl|g^<c|l8(v(j)/8x||o,/^ = c||r(o||o,/^<c‘l|rl|g^ . 

and, by the Poincare inequality, ||v|lo^ < c llrllg^^, . Using the following inverse inequality [CQ1 , 
Lemma 2.1], valid for any real numbers r and s, 0 < r < s, 

(IV.29) V4)€P,(0), , 

we can control the norm IMIg,^ of 3v/8y by 
llav/8yilo_^<cN^l|r|lo^ • 

Working out in a symmetric way for the component s of q, we end up with an element v = (v,w) of 
X^ satisfying 

(IV. 30) o’’ div (vw) = - q 

(sothatb,(v,q)=||q||J^)and 

(IV.31) ll»ll,,„<cN“(||r||,„*||s||„„) . 

It remains to estimate the gap (IV.2I) between A and B. To this end, set tp = X, Tq + Xg 
in A and ip = Tj;, + Pj in B. From (IV. 1) and (IV. 9), it follows that 
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(1V.32) b|g_^ = n(X?+X 2 ^/ 2 ) and ikilo,, = Ti (N^ yj + (N-. 1 y^) , 

Moreover, we have 

(To,Tjg = 0 if mis even and (Tq , = n m if mis odd , 

(T, ,T,)^=0 and (^ . = n (N. 1 ) . 

so that 

(IV. 33) 1 (ip, 4 /)J < Ti (|X, ||y,| N + |Xj|y 2 l (N+1) + |X 2 l|y 2 l (N+1) ) 

From (IV. 32) and (IV. 33) , we conclude that 

(IV. 34) V^€A,V.^eB, |((|).4')pNcN-^^^iklloJkllo,, . 

Since for N ^ 3 the polynomials Tq , , T,;, and are linearly independent, 6 (A,B) is greater 

than 0. Then, due to (IV. 22), (IV. 23) and (IV. 34), it is bounded from below independently of N. 
Hence, the proposition follows from (IV. 30), (IV.31 ) and (IV. 25). 

Proposition IV. 4 ; Let the space M 2 jj be defined by (IV. 1 8 ) 2 . There exists a constant I 2 > 0 
independent of N such that 

(IV.35) Vq€M 2 , . sup ^4^ ^& 2 N*'llqlloc. ■ 

v€X|,j 

Proof : In the space H = Pf^(-1 ,1) provided with the inner product (.,.)^ , we set now 
A = { Qn , Tn } and B = ( Tilj , }; due to the definition (IV. 1 6 ) of , A"'" is the subspace of 

Pn-i(- 1 .1) of the polynomials having zero average on (- 1 , 1 ). From Proposition IV. 2, we see that 
M 2 fj is equal to (A®^)'*' fl (B®^)"*". Thus, by Lemma IV. 4, each q in M 2 |,j can be expanded as in 
(IV. 24) with 

(IV.36) r(x,y)=0-u^)Z:l‘; a„(L|)L„(x) ■ « Pk- 2 (-' ■' > ■ 

and 

(IV.37) s(x.y)=(l-x^)L':J b„(x)L„(y) . b„ € P^_ 2 (- 1 . 1 ) 

(here, we find more appropriate to use a Legendre expansion for r and s). Following the proof of 

Proposition IV. 3, but using the differential equation satisfied by the Legendre polynomials (see 

also [B MM2, § V]), we define an element v = (v,w) in such that 

(IV. 38) divv = -q 

(so that b 2 (v,q) = ||q|lo_^j) and 

(IV.39) l|v||,„<cNni|r||„„t||s||„„) , 

In order to estimate the gap between A and B, let i|i = X, q^ + X 2 T,g and = y, T,;, + p 2 T,;j^, be 
arbitrary elements in A and B respectively. By (IV. 1 6) and (IV. 4) we have 
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SO that 

(Qjj , Ttr,)p = 0 if mis odd and (q^ , = 2/( 1 -m^) if mis even 

Thus||(p||o ^ ^ c (xj + Xg). Moreover, 

= 1-1 if mis even and (qN-''’m)p = 2 if m is odd . 

As in the proof of Proposition IV. 3, we have again (IV, 34), hence 6(A,B) is bounded from below 
independently of N. Thus the proposition is proved. 

Remark IV. 4 : Unfortunately, the constants involved in the inf-sup conditions (IV, 26) and 
(IV. 35) are not independent of N. However, we shall see that this does not corrupt the accuracy of 
the approximation of the velocity in the homogeneous case. 

IV. 4. The Inf-suD condition for the form a. 

Let us denote here by Kjjj (i = 1,2) the discrete kernels 
(IV.40) KjM = (v€Xf, ; Vv€ Mj^ ,bj(v,q) =0} , 

According to (lV.7)jand (IV, 1 8); , we actually have 
K,^ = { V € ;div (vw) = 0 in Q } , 

K 2 N = { V € X,^ ; divv = 0 in Q } , 

or equivalently 

(IV.41) Kj^ = K,nx^ (i = 1,2) , 

where continuous kernels Kj are defined in (III. 1 7) and (1 1 1.1 8). Moreover, by Propositions IV. 1 

and IV. 2, we have 

(IV. 42) dim = dim K 2 fj 

Thus, by Remark II. 1 , it is enough to check the inf-sup condition (11,19) for the form a. 

Proposition IV. 5 : There exists a constant oc > 0 independent of N such that 

(IV.43) Vu€< 2 n, sup ^Bcllull,^ . 

V € K ^ 1,(a) 

Proof ; Let u = (u,v) be an element of < 2 ^ . By (IV.41 ), there exists ip in H^q(Q) such that 
u = curl ip. Since 

ip(x,y) = -j“i v(^,y)d^ = J^, u(x,-o)d-o , 

necessarily ip belongs to P^(0) (see also [S]). Thus the element v = u"’ curl (ipco) used in the 
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proof of Proposition 1 1 1.2 to check the continuous inf-sup condition is indeed an element of 
K, n , whence the result. 

Remark IV. 5 ; The constant a is the constant a of Proposition 111.2, hence it is independent 
of N. 

IV. 5. A convergence estimate. 

We have proved above that the abstract assumptions (11.1 9), (11.21 ) and (ll.22)j (i = 1,2) 
are fulfilled with the present choice of spaces for the discrete velocity and pressure. Applying 
Corollary 11.2, we derive the following result. 

Theorem IV. 1 : For each integer N > 3, the Galerkin approximation (IV. 6) to the Stokes 
problem (III. 1)(lll.2) has a unique solution (u^, p*^) in x M,,g , where isdefined by 
( IV. 1 8) , . Moreover, the following inequality is satisfied 
(IV.44) ||u"||,„«N-'||p“||„^«c||l|t 

for a constant c> 0 independent of N. 

We can obtain a convergence estimate for the velocity using the abstract error estimate 
(11.35). To this end, note that by (IV.41 ) , the inclusion (1 1.34) holds. Thus we get 
(IV.45) ||u-u'^ll,^<c inf llu-vjl,^ . 

€ KgN 

We now recall a general result of polynomial approximation theory in the Chebyshev 
norms. This result is due to [Ma2] , except for the case m = 0 (where q( 0) is the space L^(Q) 
and TTJJ is the truncation of the Chebyshev series) for which we refer to [CQ1 ]. 

Lemma IV. 5 : For each integer m ^ 0, there exists a projection operator TT JJ’ : 

^ n h;;; o(0) such that, for 0 ^ r ^ m < s, 

(1V.46) Vq)€Hyo), l|.p-n,"’4)IU<c(r,m,s)N''-=||(p||^^ . 

Using the previous lemma one can show that divergence-free vector fields can be 
approximated by divergence-free polynomial fields with an optimal error estimate in the weighted 
Sobolev scale. 

Lemma IV. 6 : For each v in Kj, there exists an element Q^jV of X^ satisfying div (Q^v) = 0 in 
O such that, if v belongs to H®(0) for a real number s ^ 1 , 
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(IV. 47) 1|v-Qnv11i ^ ^ c . 

Proof ; A general proof, which applies to vector fields in IR^ as well , has been given in [SV], For 
convenience of the reader we give here a simpler proof, which however holds in the 2-dimensional 
case only. Each v satisfying the assumptions of the lemma can be written v = curl (p, with (p in 
n Define ip^ in H H^q(O) as ip^ = TT^ip and set Q^v = curl (p,,,. Then 

(IV. 47) is a direct consequence of (IV. 46). 

Using (IV. 47) , we derive from (IV. 45) the following convergence result. 

Theorem IV.2 : Assume that the solution (u,p) of the Stokes problem (111.1 )(1I1.2) belongs 
to [H^(Q)]^ X H^‘Vo) for a real number s$^ 1. Then the approximate velocity , as defined in 
Theorem IV. 1 , satisfies the convergence estimate 
(IV. 48) 11u-u'1Ii^,^cn'-M1u1I,^ 

for a constant c > 0 independent of N. 

IV. 6 Computation of th e pressure. 

Unfortunately, it is not possible to establish a similar optimal result for the pressure, if 
the discrete space of pressures is defined by (IV. 18), . As a matter of fact, it follows from 
Proposition IV. 1 and Corollary IV. 1 that the elements of M,,, are polynomials which vanish at the 
four corners of the domain, while the exact pressure needs not satisfy this condition. Thus, 
spectral accuracy in the pressure cannot be achieved. 

The remedy consists in choosing as discrete space of pressures another supplementary space 
to Z,|,j in the space P,,j(Q), which exhibits better approximation properties to the functions of M, 
and, at the same time, which fulfills an inf-sup condition asymptotically not worse than (IV. 26). 

From a general point of view, the new space of discrete pressures can be defined by 
introducing a new subspace 2,,^ of dimension 8 (which will be a suitable perturbation of Z,,,) and 
then setting 

(1V.49) = { q € P,(0) i V r € Z,, , (q,r)^ = 0 } . 

The space M,,, will allow the exact pressure to be approximated at a spectral rate if there 
exists a real number X, 0 < X < I , such that 
(IV.50) P[,^,„(Q) DM, c M,,, . 

([XN] is the integral part of XN). This means that the elements of Z,,., should be orthogonal to 
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P[,N,(0)nM,. 

On the other hand, as far as the inf-sup condition for b,^ is concerned, let -> M,,g 

denote the orthogonal projection onto M,,g with respect to the inner product (,,.)^^ . We make the 
assumption that there exists a constant c> 0 independent of N such that 
(IV.51) Vq€~M,, , llqllo.^<cl|n,qllo^ . 


Proposition IV, 6 ; Let be a subspace of P^CO) such that hypothesis (IV.51) holds. There 
exists a constant > 0 independent of N such that 


(1V.52) VqeM,^, sup 


v€X. 


bl(v.q) .2 m-2 

^ Pi 


M,u 


II0.U, 


Proof : Due to the definition (IV. 1 8), of , it follows that one has for all q in 
(IV. 53) V V 6 Xfj , b,(v,q) = b,(v,n,^q) , 


from which we deduce 


sup 
v€ X 


N 


b,(v,q) 

IHI1.W 


sup 
v€ X 


N 


b,(v,n,,,q) 

llvll,.. 


The result follows by using (IV.51 ). 




Let denote now the kernel defined as in (IV, 40) , with replaced by . If (IV.51) 
holds, K,,,, still satisfies (IV. 41), thanks to (IV. 53). Hence, both (IV. 42) and Proposition IV. 5 
are valid. Applying Corollary 11.2, we obtain the same result as in Theorem IV. 1 . 

Theorem IV. 3 : For each integer N ^ 3, the Galerkin approximation (IV. 6) to the Stokes 
problem (III. 1 )(l 11.2) has a unique solution (u*^, p*^) in x , where satisfies the 
hypothesis (IV.51). Moreover, the following inequality is satisfied 
(IV.S4) l|u"||,„.N-*||p"||„„<c||l||,. 

for a constant c > 0 independent of N. 


Remark IV. 6 : It follows from (IV. 53) that, if (u*^, p^ is the solution of (IV. 6) in X,^ x , 

then (u*^, TXfjP*^) is the solution of (IV. 6) in X^ x . In particular, the discrete velocity is 
independent of the choice of the space of pressures. 


We can obtain a convergence estimate for both the velocity and the pressure. 

Theorem IV. 4 : Assume that hypotheses (IV. 50) and (IV.51 ) hold and that the solution (u,p) of 
the Stokes problem (III. 1 )(lll.2) belongs to [H^(Q)]^x H^"’(0) for a real number s^ 1, Then 
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the approximate solution (u^, p*^) in x M,,g, as defined in Theorem IV, 3, satisfies the 

convergence estimates 

(IV.55) ||u-u"||,„<cN'-’||u||,„ 

(IV.55) llp-p"llo,„<oN’->(||u||,„*||pt,„) 

for a constant c> 0 independent of N. 

Proof : It is a straightforward consequence of (11.38) and (11.39), if we note that, due to 
(IV. 50), 

IIP-^N II 1 .U, ^ llP-'rt?XN]Pllo,u> < C ||p||^_i_^ . 

Pn^^in 

Now, we give an example of choice of the space M^|,j (or, equivalently, of Z.,|,j) for which 
(IV. 50) and (1V.51 ) hold. The argument is an adaptation of [BMM2, § V.3]. 

Let us recall that Z,^ is spanned by { Tq , T,,j and { Tilj , }®^. By (IV. 4), the 

polynomials s,,j = Tilj and t,,j = admit the Chebyshev expansions 
(IV. 57) s^ = = 2N ( + T ^_3 + ... + T^_( 2 j^D + ... ) 

Let us define the polynomials 

(IV. 59) = H\N <m<N ‘^m ^ "'"n- 3 ^mo ^ ’ 

where mg is the smallest integer > XN for which is 0 , and 
(IV. 6 O) = 6 „T„ = 2(N*I)(T„_2.T, T„^) , 

where ng is the smallest integer > XN for which is 0. Finally, let us set 

(1V.61) Z,,3Span[{To,T,)®2U{iN.tN)®'] 

and define by (IV. 49). 

Proposition IV. 7 : The space defined by (IV. 49) and (IV. 61) satisfy the hypotheses 
(IV. 50) and (1V.51). 

Proof ; The inclusion (IV. 50) holds by construction of 2,^, , due to the definitions (IV. 59) and 
(IV. 60) of i,g and t^j .Let us check (IV. 5 1 ). To this end, we shall exhibit the inverse mapping of the 
projection . For each element »i; in ( s^ , t^ }®^, we denote by the corresponding element in 
{ i|,j , t^ }®^ , and by ij; its projection onto the orthogonal complement to { Tq , T^ }®^ in L^(Q). It is 
easily seen that ij; has the form 
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*1/ = <J< - L iP . 

4>€{To.T, )®2 

for suitable coefficients . Thus, belongs to , according to Proposition IV, 1 . For each q in 
, let us define 


(IV. 62) -CMd^q 




02 




where the coefficients are such that i,gq belongs to , i.e. , 

(we have used the fact that, if iv, and are two different elements of ( s^^ , t^^ then n;, is 

orthogonal to v^g). Since »}; belongs to , we have b,(v, q-ij^q) 0 for all v in , or 
equivalently, reoalling the definition of M,,g , 

VreM,^, (q-x^Q. f')co = 0 ■ 

Thus q is the orthogonal projection of t^q onto , i.e., is the inverse of n^j , Finally, let us 
estimate the norm of X|,j . By (IV. 62) and (IV. 63), we have for all q in 

llx^qllo,^ <llq|lo,c. + ^ lYjlltllo.o. ■ 

'<' € ( Sfj , tfj }®^ 


^ llqllo,u> ^ ^ 


L 


4' € { S. , t. }®^ 


^o.w ) 

0,u> 


Taking into account (IV. 1), (IV. 57) and (IV. 58), one obtains 

while by (IV, 1 ) , (IV.59) and (IV.60) , one has 
l|iJlo,p^cNl"'2 and lit Jlo.p ^ c , 

It follows that, for any ij> in { s^ , t,,) }®^, Ikllou ^ H^Hou bounded independently of N, thus for 
each q in 

llx^qllo.^, <c||q|lo^ - 
which is precisely (IV. 5 1 ), 


Due to Proposition IV, 7, we oan apply our algorithm with the space defined by (IV. 49) 
and (IV. 61 ), and we do obtain the estimates (IV. 55) and (IV. 56). 
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V A collocation method for the Stokes problem. 

In this section, we will study the convergence properties of a spectral collocation 
approximation to the Stokes problem (111,1). This method can be extended to an efficient 
collocation scheme for the approximation of the full Navier-Stokes equations (see [Me]) ; the 
error analysis for this latter approximation will be performed in the following section. 

V.1. The discrete problem. 

We consider the simplest Chebyshev collocation approximation of the Stokes problem 
(lll.l). This scheme uses a single grid for both the momentum and the continuity equation, the 
grid being given by the cartesian product of the Gauss-Lobatto points in one space variable for the 
Chebyshev mesure p(0 dC. 

More precisely, for a fixed N > 3, let us set for 0 < j < N 
(V.l) Cj = cos(Jti/N) and Qj = Tt/c^N , 

(with Cq = cjg = 2 and c^ = 1 for 1 < j < N- 1 ), These are respectively the knots and the weights of 
the Gauss-Lobatto integration formula for the Chebyshev weight, exact for polynomials of degree 
< 2N-1 , i.e., 

(V.2) V(p€P2N.i(-i,i), (p(C)Q(c;)cic = Lj^,o'p(f;pej • 

It will be useful in the sequel to recall that the interior quadrature nodes are the zeros of the 
polynomial T,;, , i.e., 

(V.3) V j, 1 ^ j <N-1 , T;^(Cj) = 0 , 

We recall also that, as a consequence of (V.2), the bilinear form 

(V.4) V(«,»)€[tf(l-t,l])]L = *(9 6, , 

is indeed a scalar product on P^(-1 ,1) and the associated norm is uniformly equivalent to the 

norm IMIq^ since we have (see [CQl , §3]) 

(V.5) Vif €P,(-1,1), ||x)llo.e<(‘P.X'),.r<^ll'Pllo,8 • 

Finally, for each function tp in (f(Q), we shall denote by /^tp the unique polynomial in P^CO) 
which interpolates ip at the nodes defined in (V. 1), i.e. , such that 
(V.6) Vj,0^j<N, /Ntp(C;j) = -pCCj) • 

Next we introduce the cartesian product of the points defined in (V. 1 ) 

(V.7) Zfj = {x = (C;. ,cp.o^j,k<N} , 
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as well as the corresponding weights 

(V.8) Vx = (c;j .Ck) • ^x = ejek ■ 

The bilinear form 

(V.9) V (ip,.!;) € [C^(Q)]^ (41,11;)^^ = ■ 

is an inner product on P,^(Q). The associated norm, which we denote by , is uniformly 

equivalent to the norm 1 |.||q^ ( see [CQ1 , §3]), more precisely we have, 

(V.10) VtpeP^CQ), ll<p|lo.<^^ll>plU<2||4)|lo.^ . 

Finally, for each function f in (fiO), we shall denote by the unique polynomial in Pfj(O) 
which interpolates ip at the nodes defined in (V.7) , i.e, , 

(V.ll) VX€Z^. u/fj(p(x) = ip(x) , 

Note that we have 

(V.12) Vip€P^(0), (tp -o/^tp.>V)^.N = 0 ■ 
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Figure V.l 

The set of collocation nodes (for N = 7), 


Let us now introduce the collocation approximation to problem (111,1), first in the 
homogeneous case (1 1 1.2). We look for the approximate velocity in the space = [P^ (0)]^ 
and for the approximate pressure p*^ in P,^(0). We assume here that f belongs to [C^(0)]^, We 
consider the following problem : Find (u^.p*^) in X^ x Pm(0) such that 


(V.13) 


- V Au^(x) + (grad p'^)(x) = f(x) forx€ 
(div )(x) =0 for X € Zm . 


no 


In order to discuss the well-posedness of this problem, as well as its convergence 
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properties, we now give a variational formulation of (V. 1 3) which fits into the abstract scheme 
(11.17). Thus let us introduce the following bilinear forms a,g : x R and bj^ : 

X^ X P|^(Q)-». [R (i = 1 , 2) defined by 
(V.14) a,,j(u,v) = - V (Au.v)^^fj , 

(V.15) b,fj(v,q) = (v,grad , 

and 

(V.16) b 2 N(v,q) = - (div . 

Let us note that, if q belongs to P,g(0), grad q belongs to [Pn_i(- 1 .1 )®Pn(-1 ,D] x 
[P fj(-1 ,1 )®Pn_i(- 1 -O]’ Hence, from (V.2) we deduce that for any v= (v,w) in X,.j 
b,^(v,q) = q(x) dx (8q/ax)(x,t;|^) v(x,C,,) Q,, 

+ Lj^o Pj 1-1 (9q/9y)(Cj ,y) w(Cj ,y) g(y) dy . 

Let us integrate by parts the first term in the x-direction and the second one in the y-direction . 
Noting that div(vto) belongs to [Pn_i(-1 .1)(S>Pn(-1 J)] + [Pn(-1 .0®Pn.i(- 1 ,D] and 
recalling (V.2) we obtain 
(V.17) b,fj(v,q) = (u"’ div(vu),q)^,j . 

Using the same argument we also have 

(V.18) a,g(u,v) = V (grad u, u”' grad (vu))j^^ 

Thus the bilinear forms a,^ and bj^ (i = 1,2) are discrete approximations of the forms a and bj 
(i = 1 , 2) defined in (111.14) to (111.16). 


Proposition V. 1 ; Problem (V. 1 3) is equivalent to the following variational one : Find (u'^,p’^) 
in \ X ? ^{Cl) such that 
(V.19) 


Vv€ X 


N 


a^(u'^,v) + b,fj(v,p'^) = (t,v) 


u),N 


Vq€P^(Q), b 2 N(u^q) = 0 


Proof ; The first equation in (V. 1 9) is obtained by taking the inner product in of the first 
equation in (V. 1 3) by v(x) , for x in H 0 , and summing up over the points of (let us 
recall that, since v belongs to , it is equal to 0 at any point of Z,g of the boundary of D). 
Similarly, the second equation in (V. 1 9) is obtained by multiplying the second equation in (V. 1 3) 
by q(x) and summing up over all x in . Conversely (V. 1 3) follows from (V. 1 9) using as test 
functions the Lagrange basis in X,,, and P^(C)) associated to the set of points Z,,j . 
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In analogy to the Galerkin approximation, one expects that the spaces of trial and test 
pressures have to be restricted to proper subspaces and of Pn(0) so that the bilinear 
forms a^j and b^^ (i = 1 , 2) satisfy the inf-sup conditions (11.19), (1 1.20) and (1 1.22)^ . Then, we 
obtain another particular case of the abstract approximate problem (11.17), if we set, as before, 
= Xgg □ provided with the norm of X, and M 28 = provided with the norm 

||.||q^ , and if the forms a^ and bjg (i = 1 , 2) are respectively the forms a^ and b|^^ (i = 1 , 2) defined 
in (V.14),(V.15) and (V.16). 

Hereafter we shall characteri 2 e the spurious modes for the pressure and we shall indicate 
again one choice of the spaces and M 2 ,,j which leads to a "spectral" rate of convergence for both 
the velocity field and the pressure. 

V.2. The spurious modes of the pressure. 

The characterization of the parasitic modes for the discrete pressure can be carried out by 
suitably modifying the proofs given in Section IV. 2 in the case of the Galerkin approximation. 
However, we prefer to follow a different strategy, which consists in transferring the results 
therein obtained into the context of a collocation approximation. To this end, let us define the 
operators : P^(-l ,1) -► Pn^~' .1 ) by 
(V.20) V € Pfj(-1 ,1 ), (Rn‘P.'I»)p,n = (>P.4')g . 
then S|,j : P^(Q) -> Pn( 0) by 
(V.21) V4iePfj(0), (S,g(p,4<)^^N = (ip,4')ij • 

By (V.5) and (V. 10), both and S,g are isomorphisms, the norms of which can be bounded 
independently of N if we endow P^(-l ,1 ) and P^(Q) respectively with the norms |1.||q^ and |1.||q ^ . 

Now, let Zj^ denote the subspace of all q in P,,j(0) for which we have 
(V.22)j Vv€X^, bj^(v,q) = 0 . 

In order to characterize the spaces Zjjg , we introduce the polynomials Cq and r^ of Ph,(- 1 , 1 ) 
satisfying 

(V.23) V j,O^UN , ro(Cj) =60. and r^(Cj) = 6^^ , 

and a polynomial 
(V.24) q|,j* □ RfjRfg I 

where q^^ is defined in (IV. 1 6). It is an easy matter to check that q^j* satisfies the relation 
(V.25) Vip€Pn(-1,D. (qN*.4>)p,N = l^ 'l>(C)dC . 
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In analogy to Proposition IV. 1 , we can state the following result. 

Proposition V.2 : The subset Z,,^ is the vector space of dimension 8 spanned by { Tq , }®^ 

and by [fq, 

Proof : Let q in P,^ satisfy (V.22)^ . Then is such that 
V V € X|^ , b^(v.S,;j' q) = 0 , 

hence, by proposition IV.l, q belongs to the space spanned by { R^Tq . }®^ and 

{ R^Tilj , R^Tiij^, }®^. By (V.2) we have R^Tq = Tq and R,.jT|;j = T,;j , while a direct computation shows 
that (see e.g. [CQl , §3]) 

(V.26) (Tn .V^,.2(T, ,T,)^ , 

so that R^T|^ = 2 Tfj . Let us check that R|gT|!j^i(C) = ((N + 1)/N) C T||j(0- Differentiating (IV. 3) and 
using (IV. 4) yields 

(V.27) T;j^,(C;) = ((N+l)/N)(CT'(t;) +NT^(0) . 

Thus the relation 

((N+1)/N)(CT,,vi.)^, = (T,^, , 

holds for all »i( in Pfj.i( - 1 , 1 ) . as a consequence of (V.2) and (IV. 1 ). For t|; = T^ , the identity 
follows from (IV. 4) , (V.2) , (V.26) and (V.27). 

Hence we have proven that Z,,^ = Span [{ Tq , T^ }®^ U { T,;j , CT,;j }®^]. Finally we easily 
derive from (IV. 2) and (V.3) that 

(V.28) ro(C) = (-1)V(UC)/2 N^)t;j(C) and r^(C) = (( 1 -C)/2N^) T;^(C) , 

which ends the proof. 

Following the same lines we can prove the analogous result to Proposition IV. 2. 


Proposition V.5 : The subspace is the vector space of dimension 8 spanned by 
{Qn* -T^)®^ancfby {ro ,rj®^. 
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i = 1 

i = 2 

Form bj 



Form bjfj 

(T„,T„ 

{q„M„ )®'U{r„,r„ 


Figure V.2 

The parasitic modes for the forms bj and b||^ (i = 1 , 2). 


V.5, An inf-suD condition for the forms b .„ (i = 1. 2). 

In order to satisfy the discrete inf-sup condition (ll.22)| for the forms b^^ (i=1 , 2). 
according to Propositions V.2 and V.3, we oan ohoose as spaces (i = 1 , 2) any supplementary 
space to Z|fj in P^(0), i.e. any subspace of Pfj(O) such that 
(V.29) [ codim = 8 , 


1 M*nz* = {o) , 

If this condition is fulfilled, the existenoe of a oonstant §j|,j for which (II.22)j holds is ensured by 
the finite dimension of the spaces Mj|g . As in the previous case of approximation by Galerkin 
spectral method, the choice 


(V.30), h., = {q € P,(Q) ; V r € , (q,r)^, = 0 } , 

leads to a minimal constant and is useful to prepare a better choice, Alternative characterizations 
for (V.30)| follows from Propositions V.2 and V.3 : we have as in the Galerkin approximation (see 


(IV. 19) and (IV. 20)) 


(V.31) 

and 

(V.32) 


= { CO’’ div (vcj), V € \ ] 

= (qePN(0) ; V r e { Tq , (q,r)^ = 0 and q(± 1 ,± 1 ) = 0 } , 

= { q € P|g(Q) ; V r e { q,,j , (q,r)^ = 0 and q(± 1 ,± 1 ) = 0 } . 


We are now able to precise the constant for this ohoice of spaces (i = 1 , 2). 
Proposition V.4 : Let the space be defined by (V.30)j . There exists a constant 0 
independent of N such that 


(V.33), VqeM,^, sup 

V € X, 


bjii(v,q) 


> liN’^llql 




O.u) 
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Proof : Consider the case i = 1 first. We know from the proof of Proposition IV. 3 that, for each q 
in , there exists v in such that - co’’ div(vo)) = qandilvH, < c (IV. 30) and 

(IV.31 )). Since the discrete inner product induces a norm uniformly equivalent to the norm ||.||q^ 
on P^Cfi) (see (V. 1 0)) , we have 

b,N(v.q) = (q,q)^,N^c||q||o^ , 

whence (V.33), holds. The case i □ 2 follows similarly using now the proof of Proposition IV. 4 . 


As established in Section IV. 6. the choice (V.30)j is not well suited to approximate the 
space of pressures. Here again an alternative choice for the pressure space can be given. In 
order to ensure the exact pressure in M, to be approximated within spectral accuracy by an 
element in , we require that there exists a real number X, 0 < X < 1 , such that 
(V.34) . 

Furthermore, in order to retain the compatibility between and , we also require that the 
orthogonal projection Ttjjj -. ^ respect to the inner product (•,.)^iq . 

satisfies 

(V.35) Vq€M,, , llq|lo.,.<c!lTT*qllo^ . 


This yields, as in Section IV , the following inf-sup condition : 

Proposition V.5 : Let h be a subspace of ? such that hypothesis holds. There 

exists a constant 0 independent of N such that 

(V.36) Vq€M,^ 


sup 

v€X 




N 


h.u> 


IIO.u 


We refer to Section IV to the proof of Proposition V.5 and for the discussion of the existence 
of spaces satisfying the hypotheses (V.34) and (V.35); indeed the space defined by 
(IV. 49) and (IV.61 ) works (see also Section V.6 for some considerations on the implementation of 
the method). 


V.4.The inf-suD condition for the form_a„ _. 

Let us assume here that (i = 1 , 2) is any supplementary space in P^(f^) of the space Zj,^ 
defined in (V.22)j , i.e. satisfies (V.29). Setting 
(V.37) KiN = {v€XN;VqeMi^.biJv.q) = 0) . 
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by the definition of Zj^ . we actually have for all v in 
Vq€P|^(0), bj,^(v,q) = 0 , 

thus, as in the Galerkin method, 

(V.38) (i=1.2) . 

It is readily seen that in this case we have also dim = dim Kg,,, . Thus, by Remark II. 1 , it is 
enough to check the inf-sup condition (II, 1 9) for the form a^ over x . 

Prnnosition V.6 : There exists a constant a > 0 independent of N such that 

(V.39) VueKgfj, sup Yii’ ~ 

V € ll''lli,w 

Proof : Let u = (u,v) be an element of . As in the proof of Proposition IV. 5, we deduce from 
(V.38) that there exists ip in Pfj(O) Pi q (Q) such that 

(V.40) u = curl ip . 

Let us set again v r co"' curl (tpco) and recall (see the proof of Proposition 1 1 1.2) that 
(V.41) IMI,_^^c|lul|,_^ . 

Setting now, as in the continuous case, x = Aip , we have by the definition (V. 1 4) of the form a^^ 
a^(u,v) = - V (Au,v)^^fj = - V (curl x, curl (ipu))^fj 
or more explicitely 

(V.42) a^(u,v) = - V (0x/ax, co"' a(ipcj)/3x)^_^ - (dx/dy, co'' a(ipcj)/ay)^,,j . 

On the other hand, recalling that the Chebyshev weight cj is the product of the Chebyshev weights 
in each direction, i.e. cj(x) = Q(x)g(y), we have 

- (ax/ax, a(ipcj)/ax)^,g 

= -LLo (9x/9x)(Cj,c;k)[(9'()/ax)(t;j,c;k) + i>(Cj.t;k)e’(Cj)/g(c;j)]Qj} . 

By (V.2) and(V.3) it follows that 

- (ax/ax, a(tpco)/ax)^^,g =-Lk=oPktf^ (ax/9x)(x,i;|^) (a(ipQ)/ax)(x,Ck) dx] , 

= LkzoPk[l-i (9^x/9x^)(x,Ck) ip(x,Ck) Q(x) dx] , 

= Lj^o^k=o (9^X/9x^)(Cj.Ck) >p(Cj,Ck) QjQk] • 

The term in (V.42) containing the y-derivative can be handled similarly. Thus, we have proved 
that 

(V.43) a,g(u,v) = (A%,ip)^fj . 

By [MM, Lemma 3.2], there exists a constant c > 0 independent of N such that 
Vtp€P,(Q)nH2_o(Q), (A%,^)^_, >cbg^ . 



-43- 


Then, the proposition follows from (V.40), (V.41 ) and (V.43). 

V.5. A converoence estimate in the homogeneous case. 

In this section we consider the collocation approximation (V, 1 3) to the homogeneous Stokes 
problem (III. 1 )(lll.2). Let us recall that \ = [P^(0)]^. 

For an appropriate choice of the spaces of pressures and of test functions, we know from 
Propositions V.6, V.4 and V.5 that the bilinear forms a,,j , b^^g and bjfg restricted to these finite 
dimensional subspaces satisfy the inf-sup conditions (11.19), (11.21) and (1 1.22)^ (i = 1,2) for 
suitable constants. Moreover, it follows from (V. 16), (V. 1 7), (V. 1 8) and (V. 1 0) that a^ and bj^ 
(i = 1,2) are uniformly continuous over x X,g and \ x Mj^ respectively. 

Corollary 11.2 yields the following stability result. 

Theorem V. 1 : For each integer , the collocation approximation (V. 1 3) to the Stokes 

problem (111. 1 )(II1.2) has a unique solution (u*^,p^) in , where M,,g is defined by 

(V.30), (resp. in X,g x , where M,,g satisfies the hypothesis (V.35);. Moreover, the 

following inequality is satisfied 

(V.44) ||u + N Up IIo_(j < c ||f , 

for a constant c> 0 independent of N. 

Remark V. 1 : In this case as for the Galerkin method (see Remark IV. 6), the discrete velocity is 
independent of the space of pressures. 

Let us now consider the convergence of the approximation. 

Theorem V.2 : Assume that the solution (u,p) of the Stokes problem (111. 1 )(lll.2) belongs to 
[H^(Q)]^ X H^'Vo) for a real number s ^ 1 , and the data f belong to [H^(Q)]^ for a real 
number a > 1 .Then the approximate velocity as defined in Theorem V.l , satisfies the 
convergence estimate 

(V.4S) |u-u" II, „ « c ( n'-‘ ||u||,„ * N- IK II, „ ) 

for a constant c independent of N. 

Proof : Let us first remark that 
(V.46) Vve [p;.g(Q)]^ Vz€X^ , 


(a-a,g)(v,z) = 0 
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indeed, the product vz is an element of P 2 ^_i( 0 ) and the discrete integration formula in the 
definition (V.14) of a,g is exact (see (V.2)). By (V.38), we can apply Corollary 11.3 with 
Vg = Wg = Q(j_iU (the divergence-free polynomial approximation to u, the existence of which is 
guaranteed by Lemma IV. 6) to get the following error estimate for the velocity ; 


Ni 


u-u- 


ll,U) 


<c(||u-Q^_iu||i^+ sup 

2€ Xk 


II II / 


(V.47) 

Note that the term a-a,g has disappeared in this estimate due to (V.46). 

Due to Lemma IV. 6, it is sufficient to bound the error on the dataf. We have for any z in 

I - ('■ I < ic. ri. - (nS.,1. z)„ I + 1 (n2.,t, , I , 

where is the orthogonal projection operator onto [Pfj_i(0)]^ with respect to (.,.)^^ , and 
is the interpolation operator at the collocation nodes defined in (V.7). Hence we get by (V. 1 0) 


N-l'I'O.w II* ^N*llo,u ^ 11^ llo.cj 


(V.48) Vz€Xg|, 1 I < ( 

The first term on the right-hand side can be estimated by (I V.46), while the interpolation 
operator satisfies the following inequality ([CQl , Thm 3.1]), valid for any real numbers r and s, 
s > 1 and 0 ^ r ^ s. 


(V.49) V4)€H^(Q), Ik-^NUll 

This ends the proof of the theorem. 


< c N 


2r-s 


Remark V.2 : Estimate (V.45) is optimal with respect to the regularity of the solution and of the 
data. 


Let us turn now to the error estimate for the pressure. 

Theorem V.3 : Assume that hypotheses (V.34) and (V.35) hold and that the solution (u,p) of 
the Stokes problem {\\\.])i\\\. 2) belongs to [H^(Q)j^ x H^"’(0) for a real number s> 1 ,and 
the data f belong to [H“(Q)j^ for a real number a > ].Then the approximate pressure p*^ in 
, as defined in Theorem V. 1 , satisfies the convergence estimate 
(V.SO) l|p-p"||„„ <c { n’- ( llull,,, . ||p||._,„ ) ♦ Im|,„ ) 

for a constant c independent of N. 


Proof : We use Theorem 11.3 and, thanks to (V.46), we have 
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|P“P llo.w ^ 


inf (lIP-PNio.u, 

Pn ^ f^lN 


+ sup 
z€ X 


N 


+ sup 
Z€X. 


Pn^ 

“nnc 

ll-lll. 


) 


Due to Lemma IV. 6 and to (V.48), we only have to estimate the terms on the right-hand side 
concerning the pressure. Let us recall that the space we have chosen satisfies (V.34) for a 
fixed X < 1 . Thus, we can take q,g = Tt°xN]P’ truncation error is again estimated by (IV. 46). 
Moreover , using the exactness (V.2) of the quadrature formula as described in the previous proof, 
we have 

VzeX^, (b,-b,^)(z,n[°^,g]P) = 0 , 

hence the result. 


Remark V.5 ; Let us for a while consider the problem ; Find (u'^,p'^) in X^j x (resp. in 
Xn X such that 


(V.51) 


V V € X. 


3m(u^,v) + b,.(v,p^) = (f,v)^ 


'H ’ ^ “IN' 

VqeP|,j(Q), b2N(Ufj,q) = 0 

(which is problem (V. 1 9) with (f .v)^^ ^ replaced by (f ,v)^^ or equivalently with S^f replaced by f). 
Then, Theorems V. 1 to V.3 are still valid. Furthermore, since the second term in the right-hand 
side of (V.47) disappears, the estimates (V.45) and (V.50) can be replaced respectively by 
(V.52) ||u-u^| 

(V.53) |1p-p^' 

This will be used in Section VI. 


,Io^<cn'-*(||u||3^. 




) 


V.6. Concluding remarks in the h omogeneous case. 

We complete Section V.5 with some considerations on the practical implementation of the 
Chebyshev collocation scheme (V.l 3). As a matter of fact, we have indicated a pair of spaces M,^ 
and MgN fon which the collocation scheme is well-posed and guarantees spectral accuracy. Now, it 
remains to exhibit a precise set of algebraic equations, as well as of unknowns for the pressure, 
which correspond to the scheme and which are efficiently implementable. 

To this end, let S^, denote the set of the four corners of the square Q , and let S be a set of four 
collocation points in \ satisfying the following property : 

(V.54) det ( qL(xj) ) 0 , 1<J,L<4 , 
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where Xj runs through S and Ql runs through { q^j*, T|g (the polynomial q,g* is defined in 
(V.25)). 


Proposition '^.7 : Assume that hypothesis (V.54) holds. To apply the collocation scheme 
(V.13) is equivalent to solve the following linear system : Find (u'^.p^) in x .where 
is defined by (V.30), (or in X,^ x , where satisfies the hypothesis (V.35)^ such that 
(V.55) I - V Au^x) + (grad p^)(x) =f(x) for x € H 0 , 

(div u^)(x) =0 for X € Zo \ ( S U M 


Proof : Of course, any solution of (V. 13) is a solution of (V.55). Conversely, let (u'^,p^) satisfy 
(V.55). Since u'^ vanishes on 90, div is equal to 0 at the four corners of 0. Moreover, by 
Proposition V.3 , we know that p 2 fj(u^,q) is equal to 0 for all q in { q^*, }®^. By (V.55) this 

relation becomes 

Vq€{q,M,}®', (divu^(x)q(x)co, = 0 . 

Thanks to (V.54) , we obtain (div u^)(x) = 0 for any x € §. We conclude that is divergence-free 
in 0, hence (u^, p^) satisfies (V. 1 3) (and (V. 1 9)). 


As far as the choice of degrees of freedom for the pressure is concerned, it seems 
unpractical to find out a subset of oollocation points in the domain, which uniquely determine the 
polynomials of M,^ (i.e., which form a unisolvent set for M^n^' convenient to 

determine the discrete pressure through the complete set of collocation points in the domain (or, 
equivalently, to retain all the modes for the pressure). This means that the algebraic system to be 
aotually solved is underspecified. Once a solution of this system is obtained in some way, it will 
yield a "good" velocity field and a "good" pressure gradient at the collocation points (and only 
there !). In order to get a "good" pressure, i.e. the pressure satisfying an estimate like (V.50) , one 
has to extraot from the computed pressure its component along . This can be done by taking 
the orthogonal projection of the computed pressure onto with respect to (.,.)^j ■ Other 
techniques of filtering the spurious modes have been successfully applied (see, e.g. [Me]). 

V.7. The non homogeneous case. 

Let us now consider the approximation of the non homogeneous Stokes problem 
(III. 1 )(lll.39) by a collocation method. We shall suppose that this problem is well-posed, i.e. 
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that (111.40) and (111.41) hold. 

Hereafter, we assume that the space Is defined by (V.30)j and that the space 
satisfies (V.35). After the analysis of the homogeneous case, we propose the following formulation 
of the discrete problem :Find (u'^.p'^) in [Pfj(O)]^ x (resp. in [P^CO)]^ x such that 
(V.56) V V € X,g , a^(u'^,v) + b,fj(v,p'^) n (f.v)^^,,, , 

Vq€M2N. b2N(u^q) = 0 , 

u^(x)=(pj(x) for X € Sfj n Tj , J 6 2/4Z 

Remark V.4 : This formulation is not as direct as the formulation (V. 1 3) of the homogeneous 
case. Indeed, we have in mind to discretize the equations in a collocation way, hence we would like 
to satisfy each of the equations at a suitable set of points of . More precisely we would like to 

solve the following problem ; Find (u*^,p^) in [Pfj(O)]^ x (resp. in [Pfj(O)]^ x such 

that 

(V.57) - V Au*^(x) + (grad p*^)(x) = l(x) forx€E,^no , 

(div )(x) =0 for X € , 

u^x) = ^pJ(x) forxCEfjHTj, d^TU^TL . 

The last problem is clearly equivalent to the variational formulation ; Find (u^,p*^) in 
[P^(Q)]^ X (resp. in [Pfj(O)]^ x such that 
(V.58) V V € Xfj . afj(u^,v) + b,^(v,p^) = (f,v)^_^ , 

Vq€P^(Q), b2N(u^q)=0 . 

U*^(x) = ipj(x) for X € Efj n Tj , J € Z/4Z 

However , it follows from Section V.3 that the relation b 2 f^(u'^,q) = 0 can only be satisfied for all q 
in M 2 |,j and generally not for all q in Pfj(Q). Indeed, since is a polynomial, we would derive 
from (V.58) ; div = 0 exactly. In particular this imposes five conditions for at the 
boundary : 

(V.59) (div u*^)(x) = 0 for x corner of 0 , 

(V.60) Z.JS 2/42 Ifj u . rij do = 0 

These equations solely depend upon the values of at the boundary, hence upon /fjipj • In general 
these relations are not satisfied. 

Example ; Assume that m is even, 2 < m < N, and let us choose as boundary conditions 
(V.61) VJ€Z/4Z, ipj = ( L^(t;),L^(c;)) 

(note that the Legendre polynomial of degree m has by definition a zero-average). These conditions 
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fulfill the hypotheses (111,40) and (111.41). Thus, iff is chosen in X', the continuous problem is 
well-posed. But, condition (V.59) does not hold since coincide with (pj for any J € Z/4Z and 
div at a, for instance is equal to L^(-l ) - L|!^( 1 ) = - m(m + 1 ). In this example, the continuous 
problem has a solution in the weak sense only, since the divergence of the velocity is not zero at 
the four oorners. However, even if the exact velocity field is divergence free at the four corners, 
condition (V.59) need not be satisfied : in fact, the Lagrange interpolation operator will not 
generally preserve the boundary values of the first derivative. 

Hence, problem (V.57) is generally unsolvable and the formulation (V.56) is used instead. 
Besides, we can check the following proposition. 

Proposition V.8 : Anu solution of problem (V.56) in [P x (resp. in 

[Pfj(O) X M satisfies the collocation equation 

(V.62) - V Au*^(x) + (grad p*^)(x) = f(x) for x€E^f] Cl 

Remark V.5 : Due to (V.32) , solving the equation 
V q 6 M 2 N . b 2 N(u^q) = 0 

in (V.56) is equivalent to the minimization of l|div ^ condition is implemented in 

practice. We refer to [Me] for details and numerical results in the non homogeneous case. 

As in the proof of Theorem 111.2 we shall need an element in the space of trial functions (i.e, 
[P|,j(0)] ) that satisfies, in a discrete sense, the boundary condition (111,39). To this purpose we 
recall the following approximation result that can be found in [BM], 

Lemma V. 1 : There exists an operator from H^(Q) into P|g(0) such that, for any function 
ip in H^(Q) 

(V.63) vx€Z^neo, n^^(p(x) = (p(x) . 

(V.64) ll<C-”r(^ ij(p||,_^ < c ( inf II 111, w 6 2/42 ll‘P|rj~^N‘P|rj ^ 

First, we check that the discrete problem (V.56) is well-posed. 

Theorem V.4 : For each integer H ^ Z , the collocation approximation (V.56) to the Stokes 
problem (111. 1 )(lll.39) has a unique solution (u^,p^) in x M^^ , where M,^ is defined by 
(V.30), (resp. in X^ x , where M^^j satisfies the hypothesis (V,35)^. Moreover, the 
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following inequality is satisfied 

(V.65) ||u Hi + N Up IIo,oj ^ II[c*(o )]2 n II'Pj1I3/4,q ^ > 

for a constant c> 0 independent of N. 

Proof : If we define - TTi^ in , then the pair (u^,p^) is a solution of (V.56) if and 

only if the pair (u’^.p'^) satisfies : 

(V.66) Vv€X^, a,^(u^,v) + b^^(v,p^) □ - afj(TT,^^u,v) , 

VqeMgN. b2N(u^q) = - b2|g(n^,,u,q) . 

We derive the result from Theorem V. 1 and (V.64), together with the following estimate for the 
interpolation error ([CQ1 , Thm 3.1]), valid for any real numbers r and s, s > 1 /2 and 0 < r < s, 
(V.67) V ip € H^C-I.I), . 

In order to investigate the convergence properties, we shall introduce a slightly different 
approximation that satisfies condition (V.60). We first define the constant 
(V.68) ~ ^ ^j€Z/4Z Itj "Pj • 

Due to (111.41 ) and (V.67), if ipj , J 6 Z/4Z, belongs to Hp(Tj) for i > 1/2, we have 
(V.69) |e^ I < c N 2Ljg2/42 II'PjIIt.q 

Let us denote now by the vector (e^ , e^j) and remark that now /^u-e^j is a boundary 
condition that satisfies (111.40) and (111.41 ). Hence we can define the two auxiliary problems 

a) Find (u.p) in [H^(0)]^ x M, such that 

Vv€X, a(u,v) + b,(v,p) = (f,v)^^ , 

(V.70) V qe M 2 , b 2 (G,q) =0 , 

u = TTN b«J - e,g = Jn“ "®N i 

b) Find (u^,^*^) in [P,g(0)]^ x (resp. in [Pn(0)]^ x such that 

Vv€Xfj, a^(u'^,v) + b,|,,(v,p^ = (f,v)^^j^ , 

(V.71) Vq€M2N. b2N(u^q) = 0 , 

u'^(x) = ■ ®N - “ ®N ^’orx^ZfjDdO 

The error bound between the solutions of problem (III. 1 )(lll.39) and (V.56) will be obtained by 
studying the differences between u and u, u and u"^, and . 

Proposition V.9 : Assume that the boundary data (Pj , J e TUATL, belong to H^ffj) for a real 
number x ^ 3/4. The following estimate is satisfied 
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^J€Z/42 Il'f’jllt.Q 
for a constant c > 0 independent of N. 


(V.72) ||u-G11i^^cN 


Proof : It is a consequence of Theorem 1 1 1. 2 and of the approximation estimates (V.67) and 
(V.69). 


Similarly we can obtain an error bound between and . Indeed, note that is 
constant, equal to e,g . Hence, by(V.69) , we obtain the following result. 

Proposition V. 1 0 : Assume that the boundary data (pj , J e TUATL, belong to Hp(fj) for a 
real number x >3/4. The following estimate is satisfied 

(V.73) ||u^-u'^|li_^^cN'" 2 Lj€Z/ 42 li'P 

for a constant c> 0 independent of N. 


J 


In order to get now an error bound between u andu^ we first note that problem (V.71 ) is a 
discrete approximation of problem (V.70) . Hence the abstract results of Section 1 1 can be applied. 
We have the following error estimate. 

Proposition V.1 1 -.Assume that the solution (u,p) of the Stokes problem (1 1 1. 1 )(l 1 1.39) 
belongs to [H^(Q)]^ x H®"’(Q)for a real number s > 1 , that the data f belong to [H^(Q)]^for a 
real number a > 1 and that the boundary data (pj, J € TL/ATL, belong to H^(rj)for a real 


number i > 3/4. The following estimate is satisfied 

'Nil ^ ^ / ki3-s II..II . ki-tr ii,.ii ki7/2-T 


(V.74) l|S-S"ll,„<c(N^-||u||.„*N-’||t|l.,„*N"^-'Lj,„,j Ikjlll., ) 

for a constant c> 0 independent of N. 


Proof ; Let us setu = u - + 0 ,,^ . Since is constant, (u,p) is the solution in X x M, of 


(V.75) 


V V € X, a(u, v) + b,(v,p) = (f,v)^ - a(TT^|jU,v) 
Vq€M2, b2(u,q) = - b2(TTN.bU.q) • 


Next, - TT ^ (jU + e,^ is the solution of (V.66). Let us define the forms g and g^^ by the 

relations 

V q € M 2 . <g,q> = - b 2 (n^,jU,q) and V q € M 2 N , <9^ . Q> = - t> 2 N(nN tU,q) . 

Since we already noticed that the hypotheses of Corollary 11.3 are fulfilled, we derive from 
(11.35) 
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(V.76) 


||u-G^||,,,<c[ inf 

w € 


lu-wiii + 


inf 


u-v 


ve[PN.,(0)f 


1,u> 


+ sup 
Z€X, 


(f,z), - (f,z) 




M.w 


+ sup a(n^(^u,z) - a^(n^(^u ,z) 
Z€X, 


M.u 


• i,w ^ 

In opposition to what we were able to do in the homogeneous case, it seems that the only way to 


bound the term inf l|u-wH, is to use Proposition 11.1 . By Proposition V.4 we have 

w € Kg^CgN) 


(V.77) 


inf 

w € Kg^jCg^,) 




inf 


u-v 


v€ [Pn.i(0)]' 


l.U) 


g € M 2 (j 

1 ) From the definitions of g and g^ and (V.2) we write for any q in M 2 N 


sup 11 ^ 11 ] 


0,u> 


<g”g(j ) g> - ^2N^^N,b^’9^ ’ 

= t» 2 NCnfj|,u - q) - b 2 (rt^^u -n^_, i^u.q) 

Using now the uniform continuity of p 2 and b 2 ^ we derive 

I <g~gN > g> I ^ II ^N,b*^ ~^N-1.b^lll,u II 9 llo,u> ■ 


hence, by (V.64) we obtain 


(V.78) sup '9>J, 

q e M 2 N II 9 llo.u 


2) On the other hand, we compute 


inf llu-v||,.^ 

< l|u|ll.u,= I|U- 'TfN.b“ ■"®nIIi 

V € [Pn_i(0)]^ 

< llu-u||,<^ + ||u-n^b“lli.u, 


II <p 


J IIt.Q 


Thus, by (V.72), (V.64) and (V.69) we obtain 


(V.79) tnf llii-*ll,„<c(N'-||ul|.„.N’'^-'E„„, 5 lkJ|„) . 

v«lP;.,(0)l' 

3) Due to (V.48), (IV. 46) and (V.49), we have 

(V.80) sup < cN-^||f||,^ . 

z€X|,l II* 1*1, w 

4) Let us now estimate the last term in (V.76) . We have by definition of the forms a, a^j and by 
(V.2) 

1 a(n^t.^,z) - a/n^^u,z) I = I a(n^t.u - n^_, (,u, z) - a^(n^b^ - n^_, (,u, z) 1 

^ c||n,bU-n,_,,u||,^llz||,^ . 


Thus 
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< c ( n’ ® Hull 


(V.81) sup 

z€ 

Finally, estimate (V.74) follows from (V.76) to (V.81 ). 


M.u 


+ N ^J€Z/4Z ll'PjIkg ^ 


Using now the three previous propositions, we obtain the final error estimate 
Theorem V.5 : Assume that the solution (u ,p) of the Stokes prob lem ( 1 1 1 . 1 ) ( 1 1 1 , 39) belongs to 
[H^(0)]^ X H^"'(0) for a real number s ^ 1 , that the data f belong to [H^(Q)]^ for a real 
number a > 1 and that the boundary data tpj, J € TL/ATL, belong to Hp(fj) for a real number 
X ^ 3/4. Then the approximate velocity , as defined in Theorem V.4 , satisfies the convergence 
estimate 


(V.82) llu-u^ll,^ < c ( llull^^^ + N'" llfll^^ + 
for a constant c independent of N. 



We conclude with an estimate for the pressure. 

Theorem V.6 : Assume that hypotheses (V.34) and (V.35) hold and that the solution (u,p) of 
the Stokes problem (III. 1 )(lll.39) belongs to [H^(0)]^ x H^'Vo) for a real number s^ 1 , 
that the data 1 belong to [H“(0)]^ for a real number a > 1 and that the boundary data qjj, 
J € Z/4Z, belong to j) for a real number x > S/A.Then the approximate pressure in 

, asdefined in Theorem V.4, satisfies the convergence estimate 


(V.83) ||p-p"llo..<c{N^-®(||u|l ^ 


for a constant c independent of N. 


S.CJ 




11 / 2 -t 


'JeZ/42 




Proof : Using Proposition V.5 we derive from (V.56) that 


{llP-qNllo.c sup 
Pn ^ ^1N Z € X| 


+ sup 

Z € > 

Still taking = TT^^jjjP , we obtain easily (V.83). 






(f.z)c.-(f.z),.N 


1 
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We are interested now in the approximation of the full Navier Stokes equations on the 
domain Q by a collocation method. From now on, we shall denote by x = (x, , X2) the generic point 
of Q, and by w^ and Wg the components of any vector w in Given a force field f in 0 and a 
viscosity v > 0, the problem is to find a velocity field u = (u, , 02 ) and a pressure p solution of 
(VI. 1) - V Au + grad p + (u.V)u =f in Q , 

div u = 0 in Q , 

such that u satisfies the following homogeneous boundary conditions 
(VI. 2) u = 0 on 00 . 

We mean by (u .V)u the sum ^((Ou/dXj). 

Let us consider the nonlinear term in equation (VI. 1). We notice that, for any function w 
such that div w = 0 in 0 , we have 
(VI. 3) = 9(WiW)/aXj . 

The two forms are equivalent for the continuous problem, but generally not for the discrete 
problems. For reasons of numerical stability as well as to reduce the computation cost, it seems 
more convenient to choose the second expression. Hence, we set 
(VI. 4) G(w) = ^ • 

In order to study the well-posedness of problem (VI. 1 )(VI.2) in the spaces X and M, , we 
first state some properties of G. 

Lemma VI. 1 : For any f in X', the mapping G is of class C°° from [h'(0)]^ into X' and from 
X into X'. Furthermore, for any vi in X, the operator DGiw) is compact from Xinto X'. 


I 

> 


Proof ; For any u and w in X we have 

V € X, 11^, ^j ^2 lo (9(UjW)/aXj) V (0 dx I = I Jft (9(vcj)/aXj) dx | 

whence, from Lemma 1 1 1. 1 , we derive 

(VI. 5 ) Vv€X, lLigi^2 Jfi (a(UiW)/aXj)vwdx|^ cLi5ij^2 IlUiWjllo.ullvlli.u,- 
We recall the imbedding of h’^^(Q) into L^(Q) (cf. [CQ3, Thm 4. 1 ] or [BM]). Moreover, using 


theCalderbn extension theorem (see [A, Thm 4.32]) together with [G, Thm 1.4. 4.2], we know that 
the mapping : (ip,i|/) ipij/ is bilinear continuous from h’(0) x H^(Q) into h’"'(Q) for any 
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s > 0. Hence we have for 0 < e ^ 1/2 

(VI. 6) ||UjWj||o^^<c|lUiWj||i/2 ^c' llUjWjll,,^ ^c" llUill, llWjlli . 

From (VI. 5) and (VI. 6), we deduoe 

(VI. 7) VveX, I Li sig 2 In ^ I ^ > 

and 

(VI. 8) Vv€X, |Li^ig 2 In ■ 

Then, it is an easy matter to derive from (VI. 4) that G is of class C°° from [h'(Q)]^ into X' and 
from X into X'. The oompactness of DG(w) from X into X', is an easy oonsequence of (VI. 6) and of 
the compactness of the imbedding ; H^’'^(Q) c; H^^^(Q) , 0 < £ < 1/2. 


Let us recall that the forms a and b, (i = 1 , 2) are defined in (1 1 1. 1 4) to (1 1 1. 1 6). From now 
on, we identify L^(Q) with its dual space and we denote by (. ,.)^ the duality pairing between X and 
X'. As in (111.13), for any f in X' , the Navier-Stokes equations (VI. 1 )(VI.2) can be written in the 
following variational form : Find (u,p) in X x such that 


(VI. 9) 


Vv€X, a(u,v) + b,(v,p) + (G(u) ,v)^^ = 0 


[ VqeMj, b2(u,q)=0 . 

Then, we derive 

Theorem VI. 1 -.For anu iin X' , problem (,y\. 9) has at least a solution (u,p)/n Xx 


Proof : From (111.9) and the imbedding of L^(Q) into L^(Q) we first deduoe that f belongs to 
[H'’(Q)]^. Using now [GR, Chapter IV, Thm 2. 1 ] we obtain that there exists at least a pair (u,p) 
in [Hq(Q)]^ X L^(0) solution of (VI. 1), where p is defined up to an additive constant. From 
Lemma VI. 1 , G(u) is an element of X‘. Let (u',p') be the solution in X x of the Stokes problem 
with data - G(u), as defined in Theorem lll.l. Then, both (u,p) and (u',p') are solutions in 
[Hq(Q)] X L (Q) of the Stokes problem with data - G(u) ; the uniqueness of the solution of the 
Stokes problem implies that u and u' coincide and that p-p' is constant. Now , fixing the oonstant 
over p so that (III. 38) is satisfied, we see that (u,p) belongs in fact to X x M, and satisfies (VI. 9). 


Now, using Theorem lll.l , we define the operator T from X' into X : f -> u = Tf, where (u,p) 
is the solution of problem (1 1 1. 1 )(l 1 1.2). Clearly, if (u,p) is solution of the Navier-Stokes 
equations (VI. 9) then u is a solution in X of 
(VI. 10) u+TG(u)=0 . 

This formulation will be very useful in the sequel, together with 
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I emma VI. 2 : For any real number q > 2, there exists a constant c{q, \>)such that , if a solution 

(u,p) of problem (VI. 9) satisfies 

(Vl.11) ||u|lLq(n) <c(q,v) , 

the operator 1 +TDGiu) is an isomorphism of X. 

Proof : Bg the compactness result of Lemma VI. 1 , the operator 1 + TDG(u) 1s an isomorphism of 
X if and only if it is injective, i.e. the only solution (w,r) of the following linearized Stokes 
problem 

Vv€X, a(w,v) + b,(v,r) + (DG(u).w,v)^ = 0 , 

VqeM2, b2(w.q)=0 , 

is (0,0). From (111.1 9) and (VI. 5) we obtain 

va||wl|,^<cLui,j^2Ni%llo.u, • 

Next, using the imbedding of h’(Q) into any L®(Q), s < +oo, and Lemma III.1 , we have 

II % llo.w ^ II llo < II IIl<i( 0 ) II Wj I1 1 ^ c(q) c' II u j IlLqj^) II Wj || , ^ , 

so that 

voc||w||i^< c'(q) ||u|Lq(( 5 ) ||w|l, , 

and the lemma is proved with c(q,v) < va/c'(q)- 

VI. 2. A collocation method for the Navier-Stokes equations. 

We are going to introduce a collocation problem to approximate the Navier-Stokes 
equations, by using the same nodes as in the linear case. Let us recall that, for a fixed integer 
N ^ 3, X^ = [Pn( 0)]^. Henceforward, we still assume that (i = 1 , 2) is defined by (V.30)j and 
that satisfies the hypothesis (V.35). 

Before writing the problem, let us consider the nonlinear terms in these equations, i.e. the 
function G defined by (VI. 4). Obviously, if a function w of class C“ is known only by its values at 
the nodes x in 3,^ (see (V.7) for the definition of Z^j), it is easy to derive the values of w^w at the 
same nodes, whence J^(WjW). The pseudo-spectral approximation a(WjW)/9Xj consists in 
differentiating this interpolation function i.e. , to compute (WjW)/9X| . 

Assume now that the force f is given in [C°(Q)]^. Due to the previous remark, the 
collocation problem we are going to analyze is the following one ; Find (u'^,p'^)/n X,^ x such 
that 
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(VI. 12) 


- V Au'^(x) + (grad p'^)(x) + (9 ('^n ( u|^u'^))/9Xj)(x) =f(x) 

for X e 

(div u*^)(x) =0 for X € Z|^ 


no 


In order to study problem (VI. 1 2), we set for any function w of [C°(0)]^ 


(VI. 13) G,(w) = Li,i^2 S^(9 (-'n( - V - 

where is defined in (V.21 ). Then, we have 

Vv^eX^, (GN(w),v^)^ = (Li^i^2 .v^)^^-(f,v^)^^ 

and due to (V.2) 


(VI. 14) Vv^eX^ 


= - (^1 «i«2 WNjWN , w’ 8(C0V^)/aXj)^^ - (f.V^)^,N . 

Hence, the next result can be proved exactly as for the linear case. 

Proposition VI. 1 : Problem (VI. 1 2) is equivalent to the following variational one ; find a pair 
(u'^.p^);n X,^ X M,|^ such that 
(VI. 15) 


V V € X 


N 


V q € M 


2N 


a^j(u^,v) + b,fj(v,p^) + (Gfj(u'^),v)^^ = 0 

b2N(u'^,q)=0 . 


Now, by Theorem V. 1 , we can define the operator from X‘ into X,g : f ^ = T^f , where 

(u*^ ,p^) is the solution of problem (V.5 1 )(recall that it is exactly problem (V. 1 9) with (f ,v,,,)^_^ 
replaced by (f .Vfj)^^). Clearly, problem (VI. 1 5) implies that is a solution in X,,^ of 
(VI. 16) u^+T^G^(u^=0 . 


We begin by stating some results about the operator . 
Proposition VI. 2 : For any i in X', the operator satisfies 


(VI. 17) ||T,f||,_^<c 


and 


sup 


'Nlll.u 


(VI. 18) ||(T-T,)f||,_^ = 0 ; 

Moreover, if the solution 7 f belongs to [H^(Q)]^ for a real number s^ 1, it satisfies the error 
estimate 


(VI. 19) ll(T-T^)f|li,^<cN’-MlTf||^^^ 


Proof : By Proposition V.6, we obtain at once 


I 
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1,u> 


< c sup 


llv. ■ 


v„€X„ ii'Nni.i^ 


sup 




''n ^ 


'Nil 1.00 


which yields (VI. 17), Next, we have already noticed (see Remark V.3) that the error estimate 
(VI. 19) can be proved exactly by the same way as (V.45) for any f such that Tf belongs to 
[H^(0)]^. Finally, (VI. 1 8) holds by classical arguments using (VI. 1 9) and the density of [D(Q)]^ 
in X. 


To study problem (VI. 16), we shall use afixed point theorem due to M. CROUZEIX (see [C, 
Th. 2.2]), which is a refined form of the discrete implicit function theorem of [BRR]. For the 
reader's convenience, let us recall this theorem : we consider a c’-mapping from a Banach 
space X^ into itself and we assume that ujj is a point in X^^ such that DF,^(ujJj) is an isomorphism of 
X,g . We denote by and /-^(T)) , t) > 0, the quantities 

(VI . 20) llx^ Kn = IKDFjj(uJJ)) llt(X|.j,XN) 

/.„(-q) z: sup ( ||DF^(w„)-DF^(u*) 11 ^x^,Xn) ; "n ^ ^n ^^d llw^-u; < T] } . 

Theorem VI. 2 : Let us assume that 2^^ ^ T5 ^ such that 

Xn ^ ' there exists a unique solution of the equation Fn,(u'^) = 0 in the ball 

Bn = { Wn € X^ ; |w„-Un |1x|^ < T) )■ This solution satisfies 
(VI.21) Vw^eB^, IIu'^-WnIL < • II'"n(“n)IL ■ 


We are going to apply this theorem to the mapping f^ = 1 + 1,^0^ . In the sequel , we always 
assume that the function f is given in [H^(Q)j^, o > 1. We consider a solution u of the 
Navier-Stokes equations (VI. 1 )(VI.2) which is nonsingular in the following sense : the operator 
1 + TDG(u) is an isomorphism of X (by virtue of Lemma VI. 2, such solutions exist for f small 
enough !). Even in the classical Sobolev spaces, regularity results of the solution (u,p) as a 
consequence of the regularity of f are not easy to derive (see [G, §7.3]), whence we shall assume 
in the sequel that there exists a real number s > 1 such that (u,p) belongs to [H^(0)]^ x H^"’(Q). 

Let us denote by N‘ the integral part of (N-D/2. We choose for ujj the image of u by the 
projection operator from [H^q(O) D H^(0)]^ onto [9^.(0)]^ (this definition of seems very 
complicated, but the fact that Un*Un belongs to [Pn_i( 0)]^ will make the estimates more 
straightforward, as it will appear later). It has been proved in [Ma2] that the following estimate 
holds for any real number r , 0 < r < s, 

(VI.22) l|u-u*||,.^^cN''-^||u||^^ 
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VI.5. Some properties of the mappinQ F „. 

In order to bound the constants , /.^ and , we need several lemmas. 

I nmma VI.5 : For any real number e > 0, there exists a constant csuch that , for any if and >i; 

in P|^(0), the following estimate is satisfied 

(VI. 23) IK 1 ^ II'PIIi.w Il’t'lli.w • 

Proof : Recalling the definition of N' and from (V. 1 2) we derive that 

so that 

+ ll(4)-n;,p)n;.i;llo^^lljj(4.-n;4')(vi^-n;4<)]llo... 

From (V. 1 0) and (V. 1 1 ) we derive 

IK 1-^n) llo.u> ^ IK»p-T(^.(p)(iif-n,J.t)llo,^ + ll'n,J.(p((j(-n^.vi;)llQ ^ 

+ ll(<p-n;4')n;4'llo.,. + IK<P-n;.p)(.i/-n^.4/)]t, 

+ lin^.tp(H(-n^.ii()ll^^jg + ll(ip-TT,J.(j))n,J.4(||^jg 

< c [||<p-n,J.(p|lLoo(f)) lk-‘rTN'*t»llo.u, 

+ I|TTn’(P|Il~(o) llo.u IIT^n’4'IIl°o(o)] 

^ [ll*P IIl°o(o) ■*■ II ^N'*P llLoo(fi)] ll*l*~^N'*t*llo,<.j 

+ ||(p-nfj.(p||o^^ ||nfj.n(||Loo(Q) 

Using now (IV. 46) and the imbedding of H^*'^^(Q) into L“(0) for any e > 0, we obtain for any 
r ^ 0 

(VI. 24) IK 1 -^N)(<Pt)|lo,w^ ^ [(ll‘f’llue/2 II 'P Hp.u II ^ N''i'llut/2] 

We derive the lemma as an easy consequence of the inverse inequality (IV. 29) and of (IV. 46). 

We can now state the following result. 

Proposition VI. 3 : For Hlarge enough , the operator DFfj(u^) = 1 + T,^DG,g(ufJ) ;s an 
isomorphism of and is bounded by a constant 'i independent of N. 

Proof : We write DFjj(uJ) in the form 

(VI. 25) DF^(u;) = [1 +TDG(u)] - (T-T^)DG(u) - T^(DG(u)-DG(u*)) 

- t^(dg-dg,)(u;;) . 


I 
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Let Wfj be any element of . Since 1 +TDG(u) is an isomorphism of X, there exists a constant Cq 
independent of N such that 

(VI. 26) 11[1 ■^TDG(u)].wJ|i_^>cJ|wJ|i_^ . 

It remains to bound the three other terms in (VI. 25). 

1 ) It follows from (VI. 1 8) and from the compactness of the operator DG(u) (see Lemma V 1.1 ) that 

IKT-T„)DO(u)||,„,| = 0 . 

Hence, for N large enough, one has 

(VI.27) ll(T-T,)DG(u).wJl,^<(Co/4)llwJli^ , 

2) It follows from (VI. 17) and from the continuity of the operator DG (see Lemma VI. 1) that 

( DG(u) — DG(u^ ^ ^ II/.(X X) ^ DG(u) — DG(U| g ) x') ^ ^ II H i w 

From the convergence of to u (see (VI. 22)), we deduce that, for N large enough, 

(VI. 28) l|T/DG(u)-DG(u*)).wJ|,^<(Co/4)||wJ|,^ . 

3) By (VI. 4), (VI. 14) and (VI. 17), we know that 

II T^(DG-DG^)(ujJ) = Hi til*' ^ tWfjjUjJ) IIq ^ + IK 1 -J^)(u^*W|^) ||q_^] 

From (VI. 23), we derive for any e > 0 

II Tfj(DG-DGfj)(ujJ).Wj^||i < C N II^n Hl.to H“nIIi,u > 
whence 

(VI. 29) II T,(DG-DG,)(u*).wJli^ < c N*'* ||u|K^ ||wjli^ . 

Finally, we conclude from (VI. 25) to (VI. 29) that, for N large enough, 

V € X,^ , llDpNtuJJ|).W|^ IK > (Cq/ 4) ||w^ 111 , 

which proves the proposition. 


Lemma VI. 4 : The constant satisfies 
(VI. 30) L,g(T))<CT) . 

Proof : Let be any element in X,g . From the linearity of DG^ and DG, we write 
||DF^(Wfj)-DF^(u^) ll;.(X|g,Xf^) = ll'*'Nt^^Nt''''N"'^N)^ll/.(X,,j,X,^) 

< ||t^ (DG(w^-uJJ))IK(X|^Xfj) II''’n llaxi^.x^) • 

Using (VI. 1 7) and the continuity of the operator DG (see Lemma VI. 1 ) yields that 
I|Tn(DG(w^-uJJ))|Kq^^_i^^) ^ c llDG(Wfj-u^) |(x,g,x') ^ ^ H''''n“‘'*n Hi,cj 
O n the other hand, as in the proof of Proposition VI. 3, we know that, for any in X|^ , 
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II T^(DG-DG^)(w^ Ufj ) . Z n II 1 

= Zli gig2 ^ [^^Ni~'-'Ni ^^N^llo.co "'■ ^ ~'^N^ ^ ^ Ho.to ^ 

SO that, by (VI. 23) and for any e > 0, 

I|Tn (DG-DG,g)(Wfj-uJJ) ll^O^i^ Xfj) ^ ^ I|v''n"UnIIi_co 
T hese two inequalities, together with the definition (VI. 20) of L^ijj) , imply (VI. 30). 

I emma VI. 5 : The constant e^satisfies 
(VI.31) £. ^c(u) +c(1) N"'" . 


Proof : From (VI. 1 0) we write F|,j(u^) in the form 
Fn(u,*)=u*.T,G,(u*)-u-TG(u) 

= (u*-u) + (T^-T)G(u) + (G(u*)-G(u)) + (G^(u*)-G(u*)) , 

which gives 

(VI. 32) ^I|Tn(G(u)-G(u*))||,^ 

.||T,(G(u;)-G,(u,*))l|,_^ 

It remains to estimate these four terms. 

1 ) Using (VI. 22) yields 

(VI, 33) I1u-u*||i_^<cN'-^1|u|1^_^ , 

2) It follows from (VI. 1 9) that 

l|(T-T^)G(u)|li^^ <cN^-M|TG(u)||3_^ . 

whence, thanks to (VI, 1 0) , 

(VI. 34) l|(T-Tf^)G(u)lli^,<cN^-"||u||3^ , 

3) Due to (VI. 1 7) and to the continuity of G (see Lemma VI. 1 ), we have 

||TjG(u)-G(u*))|li ^ c||u-u*||,_^ , 

so that 

(VI. 35) ||T^(G(u)-G(u*))|i,^ 

4) From (VI, 1 7), we derive 


< 0 ® llu| 


|T„(G(u;)-G„(u„-))||,„<C sup 


“k « Xn 


'Nlll.co 


Using the definitions (VI.4) and (VI. 1 4) of G and G,,j , we have for any in 
(G(u*)-G,(u*),v,)^ 

= - (U^X- 9(CJV^)/axp^] 

-(f.V^)^+(f.VN)^,N • 
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But, since u^*ujj belongs to [P^.iC^)]^ arid co ^ a(cov^)/9Xj belongs to [P^jCO)]^, the exact and 
discrete scalar products coincide. Hence, we obtain from (V. 48), (IV, 46) and (V. 49) 

1 (GCup-G^Cu;;) ,v^)^ 1 = I (f ,Vn)^ - (f ,Vn)^,n I < c 11 f 11^^ llv^ llo^ , 

so that 

(VI.36) 11VG(u*)-Gn(uP) 11,^< cN-'llfll^^ . 

Finally, we derive the desired estimate for e^j from inequalities (VI. 32) to (VI.36). 


VI. 4. Existence result and convergence estimates. 

We can now prove the main result of this section. 

Theorem VI. 5 : Assume that there exists a solution (u,p) of the Navier-Stokes equations 
(VI.1)(VI.2) such that the operator 1 +TDG(u) is an isomorphism of X ; assume moreover 
that it belongs to [H®(Q)]^ x H®”'(Q) for a real number s > land that the data f belong to 
[H“(Q)]^ for a real number a > 1 . For N large enough, problem (VI. 1 2) admits a solution 
(u'^,p'^) in X|^ X M,|^ where is defined by (V.30), ( resp. in X^ x , where 
satisfies the hypothesis (V.35)). Moreover, the approximate velocity satisfies the 
convergence estimate 

(VI.37) IluV 11,^ < c(u) n’-® + c(f) N-" 
for constants c(u) and c(i) independent of N, 


Proof : Using Proposition VI. 3 and Lemmas VI. 4 and VI. 5, we notice that 
2 Xn ^n^^Xn^n) < 2y . ( c(u) + c(f) N'" ) , 
so that we have : 

and the assumptions of Theorem VI. 2 are satisfied for N large enough. Hence, for each tq such that 
2 y /-nCti) < 1 . 

there exists a unique solution of (VI. 1 6) in the ball = { w,^ 6 X,g ; Hw^j -“nIIi.u, < ^ Next, 
from (VI.21 ), we derive the estimate 

11^ ^ lll.w ’ 

which, together with (VI. 22) and Lemma VI. 5, yields (VI.37). 

Next, by Proposition V.3 (resp. Proposition V.5), there exists a unique p"^ in , where 
is defined by (V.30), (resp. in Pl,,g , where Pi, satisfies the hypothesis (V.35)) such that 
Vv€X,j, b,^(v,p^) = - a,J(u^v) - (G|^(u'^),v)^ 
and the pair (u'^,p'^) is a solution of problem (VI, 12). 
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Remark VI. 1 : The error bound we obtain is the same as for the Stokes problem ; it is stil 
optimal with respect to the regularity of the solution and of the data. 


In order to state an error bound for the pressure, we need a lemma. 
i VI. 6 : Th8 approximate velocity u^, as defined in Theorem VI. 3, satisfies 


(VI. 38) sup 


(G(u)-G (u ),»,)^ ^ ^u. . 

V„ L ,, 


Proof : Letv^ be any element in . We compute 

(G(u)-G,(u^,v,)^ = (G(u)-G(u"‘),v,)^ (GCu^-G^Cu^ , v,)^ 
Lemma VI. 1 and (VI. 37) give at once 

I (G(u)-G(u^),v^)^ I < ( c(u) n’"® + c(f) N"“ ) llvfjll, 

From the definitions (VI. 4) and (VI. 14) of G and G^j , we obtain 

I(g(u")-g,(u"),»,)„i<[L,5„2 Iin-^H)(UV)II„,„1IKII,,„* 

and using (VI. 24) , (IV. 29), (V.48) , (IV. 46) and (V.49) yields for e > 0 
I (G(u^-G^(u^,V|^)^ I < ( c(u) N®'® + c N'® l|f 11^^ ) . 

Finally, these two bounds imply (VI. 38). 




Theorem VI. 4 : Assume that hypotheses (V.34) and (V.35) hold and that there exists a solution 
(u,p) of the Navier-Stokes equations (VI.1)(VI.2) such that the operator 1 + TDG(u) is an 
isomorphism of X ,• assume moreover that it belongs to [H®(Q)]^ x H®"'(0) for a real number 
s > 1 and that the data f belong to [HJO)] for a real number a > 1 . Then, the approximate 
pressure p*^ in , as defined in Theorem VI. 3, satisfies the convergence estimate 
(VI. 39) ||p-p^llo,^ ^ c(u,p) N^-® + c(f) N^-'' 
for constants c(u,p) and c(f) independent of N. 


Proof : Let us introduce the solution (u^,p^) in X,g x of the following problem : 
(VI. 40) I Vv€X^, a,(Q^v) +b,(v,p^ +(G(u),v)^ = 0 , 

Vq€P^(Q), b2(u^q)=0 . 

Since is equal to - T,gG(u) , we deduce from (VI. 1 9) that 
(VI. 41) llu-u^ll, ^ c(u) n’"® , 
moreover, we obtain (see Remark V.3) 

(VI. 42) ||p-p"|L <c(u,p)N^-® . 
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Next, due to (VI. 40) and (VI. 1 5), we notice that, for any v,g in X,g , 

b,(v^ , p^-p^ =a^(u^VN) +(G(u),v^)^ -a^(u^v^) - (Gfj(u^,v^)^ 
so that, from Proposition V.5, we deduce 

(VI.43) IIp“-p"|L„<cN^ sup 


''n ^ 


'N"l.u 


Let be any element in . By the uniform continuity of a^ , we have 

a„(u"-a“,v„) <c|lu"-i"||, J|»Jl, „<c(llu-ii"||, „ *|lu-u"|,,„) II.JI, 

SO that 


(VI. 44) ajj(u^-u^,v^) ^ ( c(u) + c(f) N ° ) ||vj|, 

Now using (VI. 44) and (VI. 38) in (VI.43) yields 
(VI. 45) ||p^-p^llo^<c(u) N^-%c(f) N^-^ , 

which, together with (VI. 42), gives (VI. 39). 
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